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ABSTRACT 

The cavitation performance of a seven-bladed, highly 
skewed propeller for a Naval auxiliary (AO-177) is 
evaluated using lifting surface numerical hydrodynamic 

methods. An important contribution to accurate predition 
of cavitation performance is shown to be an accurate model 
of the effective wake. 

A new model of the effective wake is described. A 

program to calculate unsteady time-averaged, but not 
circumferentially averaged, field point velocities is 

presented. The p r o p e 1 1 e r - i n d u c e d velocities, along with 

the original nominal wake, are combined with Huang's 
axisymmetric effective wake scheme in pie-shaped wake 
segments to determine the effective wake. 

The method is used to predict severe cavitation 

extent. This prediction is confirmed by S3PA experiments. 

Thesis Supervisor: David V. Burke 

Title: Professor of Naval Architecture 
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INTRODUCTION 



The David W. Taylor 
Development Center has 

seven-bl aded , highly skewed 
Auxiliary (AO-177). This 

concerning the state of t 
propeller design. 



Naval Sh i p 
recent ly 
propel ler 
design is 
he art in 



Research and 

designed a 

for a Na v a 1 
instructive 
compu ter-a idea 



The history of the propeller design is chronicled 
in reports by Boswell (1), Valentine and Chase (2), and 
Hendrican and Remmers (3). Computer-aided propeller 
applications included lifting line theory, used to 
determine the radial load distribution and the radial 
hydrodynamic pitch angle, and lifting surface theory, 
used to determine the final geometric pitch 
distribution, camber distribution, and final propeller 
offsets, including fillets, trailing and leading edge 
details, additional thickness added to the trailing 
edge, tip geometry, rake, and hub details. Cavitation 
performance was predicted by the method of Burrill ana 
Emerson (4) and was also evaluated by experiment. 



However, sea trials 
severe unsteady propeller 
the propeller indicated 
significantly worse than 
Burrill and Emerson method, 



of the AO-177 indicated a 
force problem. Inspection of 
that it was cavitating 
had been predicted. The 
originally designed for 
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rour-to-six-biaded 



an 



expanded 



area 



ratio o 
AO- i 77 
even 
Further 
tests. 



f 0.60, did 
propeller wit 
t hough both 
, scale e £ £ e c 



propellers with 
not extend well to the seven-b Jaded 
h an expanded area ratio of 0.77, 

were for merchant-type hulls, 
ts could have thrown off the model 
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A 



I I . 



I MP ACT OF EFFECTIVE WAKE 



The AO-177 propeller performance was evaluated at 
MIT using some numerical propeller analysis computer 
programs developed at MIT. Kerwin and Lee (5) describe 
MXT-FUF-2, a computer program to calculate forces, 
steady and unsteady, generated by a propeller. 
MXT-FPV, a steady field point velocity program, is 
described in a report by Min (6). Lee <7) reports on 
MXT-PUF-3, a computer program to predict steady and 
unsteady propeller cavitation extent and volume. All 
three programs . model the propeller as a arid of 
discrete vortex segments. The wake is also modelled as 
a vortex arid in each program. 



Cavitation performance calculations, more than 
propeller force and induced velocity calculations, are 
highly sensitive to the wake data used as input. 
Experience at MIT and elsewhere has shown that using 
the nominal wake data to approximate inflow conditions, 
wh ile acceptable for propeller force and induced 
velocity calculations, is not acceptable for cavitation 
performance calculations. The wake velocities measured 
behind a model in the absence of any propeller (the 
nominal wake) must be modified to account for change in 
boundary conditions imposed by the presence of a 
propeller. This new, modified wake is generally called 
the “effective wake. 1 ' There is, unfortunately, no 
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general method for 



anticipating the effective wa 3c e 



generated b y a 
nom i n a 1 wa k e . 



particular propeller 



i n 



aensra l 



Huang and others (8,9) have developed both theory 
and numerical schemes for calculating the effective 
wake of a propeller in an axisymmetric nominal wake. 
Det Norslce Veritas has attempted to adapt Huang's 
method to n o n - a x i s vmme z r i c wakes, such as surface ship 
wakes, by dividing the nominal wake into pie-shaped 
segments and then applying the axisymmetric wake 
calculation within each pie segment, using the nominal 
wake and the steady propeller induced velocities. This 
approach is an improvement over blindly applying the 
axisymmetric effective wake calculation to a surface 
ship's wake, but, owing to the very sharp AO-177 wake 
(even by surface ship standards), it was judged that 
some further refinement of Huang's method was 
required. 



Huang's method assumed an axisymmetric nominal 

/ 

wake and a steady propeller-induced velocity field. 
The Det Norsks Veritas method assumed a general nominal 
wake and a steady propel I er-induced velocity field. 
Short of a completely general effective wake 
calculation, the most obvious refinement was to assume 
a general nominal wake, but allow an unsteady 

propeller-induced velocity field. This was the method 



3 



used for this paper. 



Th e current method of effective wake calculations 
can be s u mm ar i zed in the following steps: 

1) Use the nominal wake and propeller data to 
calculate the unsteady propeller-induced velocity field 
at representative points both inside and outside the 
propeller disk. 

2) Divide the nominal wake into small pie-shaped 
segments as in the Det Norslce Veritas method. 

3) Calculate the effective wake velocities within 
each pie segment using the nominal wake velocities and 
the unsteady propeller-induced velocities, assuming 



that each pie 


s e gmen t 


a c t s 


as though it 


we r e 


part 


o f 


a n 


axis ymrne trie 


wa k e 


with 


the 5 ame 


r a d i a 


1 v e 


1 oc 


i ty 



distributi on . 



y - 



III. 



METHOD OF EFFECTIVE WAKE CALCULATION 



The two new aspects of the current effective wake 
calculations are the propeller-induced unsteady field point 
velocity <UFPV) calculation and the use of these 
time-averaged out not circumferentially-averaged velocities 
in a pie segment nominal wake modification calculation 
(PIEWAKE). Each of these two new methods will be 
illustrated by a typical calculation. 



The unsteady field point velocity (UFPV) program is 
basically a generalization of the steady field point 
velocity program (MIT-FPV) by Min. (6). The unsteady 
version uses the full wake lattice arrangement also 
employed in MIT-PUF-2. Program inputs are the FILE14 DATA 
file output from MIT-PUF-2 and interactive directions. 
Program outputs are either single point calculations 
directed to the user's terminal or a data file to be used 
as input into PIEWAKE, for effective wake calculations. 



UFPV first reads the FILE14 DATA input, containing 
mostly propeller geometry, singularity strengths, and 
singularity geometry. It should be noted that the number 
of propeller revolution time steps specified in MIT-PUF-2 
must be divisible by the number of propeller blades For 
effective wake calculations, the wake input to MIT-PUF-2 
ought strictly to be the effective wake; however, using the 
nominal wake to calculate the p r o p e 1 1 e r - i nd u c e d velocities 



is a good approximation. 



Th e user 
field point 
coordinates or 
grids. 



next inputs field point coordinates The 
position may be entered in cylindrical 
in terms of the rotating propeller and wake 



Th e velocity 


induced 


by one blade 


a t 


each t i me 


t hen calculated . 


At each 


time step, 


the 


vortex la 



assigned the strengths corresponding to the time 
the propeller rotation. The velocity induced 
element of the vortex grid, both blade and wake, 
s u mm ed to obtain the total velocity induced by one 
that particular angular position. 



step is 
t t i c e is 
step of 
by each 
is then 
blade in 



The velocities induced by each blade of the propeller 
are then summed to produce the total propeller-induced 
velocity for that position of the propeller. Velocities 
are calculated for different propeller positions and then 
harmonically analysed. The zeroeth harmonic, then, is the 
time-averaged p r o p e 1 1 e r - i n d u c e d velocity for that field 
point. (By averaging all points at a given radius we could 
obtain the circumferentially-averaged, time-averaged, 
propeller-induced velocity, the output of the steady 
M1T-FPV . ) 

For effective wake calculations, such velocities are 
computed for up to 60 points per radius for radii varying 
from propeller hub to 1.7 times the propeller radius. 



Between the propeller hub and propeller tip, the velocities 
are calculated at the leading edge panel on the propeller 
grid. 



For the AO-177 propeller 
radius for eight radii were 
unsteady field point velocities. 



analysis, 56 points per 
used to create a file of 



PIEVAKE is a program to perform Huang's axisymmetric 
effective wales contraction calculations in each of many 
pie-shaped segments of the wake. It uses nominal wake data 
and unsteady propeller-induced velocity data from UFPV as 
inputs to produce an effective wake velocity data file as 
output . 



After reading nominal wake data and 
velocity field data, PIEVAKE extrapolates the 
data to the hub. The innermost radius of 
velocity field is taken as the hub. (UFPV 
its innermost radius, the hub radius used in 
passed to UFPV by way of FILE14 data file.) 



the' induced 
n ora i n a I wa k e 
the i nduc ed 
had taken, as 
MIT-PUF-2 and 



In each pie segment, Huang's method is used to 
calculate an effective velocity and the effective radius 
corresponding to that effective velocity. In the numerical 
form, this method is iterative in nature. 



A finite difference equation, given in reference (9), 
to calculate the effective velocity. Given 
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i s used 



calculate the 



the 



nominal wake at radii (assumed, for the first iteration, to 
be the nominal radii) in the presence of a propeller, an 
effective velocity can be calculated if an effective 
velocity at the next radius away from the hub is known. 
The effective velocity at the outermost radius is ass ume d 
to be equal to the nominal velocity there, enabling all 
effective velocities to be calculated. 



A second finite difference equation, also given in 
reference (9), is used to calculate the effective radius 
corresponding to the set of effective velocities just 
calculated. This equation uses nominal radii, nominal 
velocities, and apparent velocities. (An apparent velocity 
is the sum of the- effective velocity and the 
propel ler- induced velocity at a point.) Given one 
effective radius, the next radius away from the hub can be 
calculated. Since the innermost radius of all velocity 
fields is set at the hub radius, ail effective radii can be 
determined sequentially from innermost to outermost. These 
radii can then be used in the finite difference equation to 
calculate the next iteration of effective velocities. 
PIEWAKE iterates five times for each pie segment. 



In addition to the effective velocity wake field, 
PIEWAKE also computes the mean velocity, the volumetric 
mean velocity of the nominal wake, the volumetric mean 
velocity of the effective wake, and the effective blockage. 



Th e effective 



velocity field then then be processed 



and used 



by MlT-PUF-3 as any other waJce field. 
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I V . 



AO-1 77 



EFFECTIVE WAKE CALCULATIONS 



Append: s I contains plots of nominal; induced, and 
effective velocities varying circumferentially at several 
specified radii. Append i s II contains the same 
information, except that the nominal and effective 
velocities are normalized on the respective volumetric mean 
velocity. These plots confirm that the modifications to 
the nominal wake outlined in Chapter III are indeed 
reasonable. Several salient conclusions can be drawn from 
these plots: 

1) Effective velocity modifications of the nominal 



v e 1 o c 


i t i es 


are 


in the expected 


direction. 


Physical ly, 


t h e 


imp o s 


i t i on 


o f 


the propeller on 


t he nomi na 1 


wake draws 


mo r e 


wa ter 


into 


the 


s 1 i p s t r earn , thus 


increasing 


the veloclt 


i e s . 



This “faster water" will come from radii greater than the 
field point being sampled; in other words, the wake will 
seem to “contract.* 1 

2) Effective velocity modifications of the nominal 
wake velocities are greatest at lower radii. The wake 
"contraction" will seem to concentrate here. 

3) In general; the wake peak is more narrow in the 
effective wake than in the nominal wake. In the case of 
the AO-177 nominal wake, the existence in the presence of a 
propeller of a wake peak as sharp as that in the nominal 
wake is counter-intuitive. Therefore, a decrease in the 
wake peak is to be expected. 

4) The normalized effective velocities at 0 degrees 



are not appreciably different from that of the nominal 
wake. This is disappointing and was not expected. It may 
be a result of radial vorticity, ignored by the 
ax i symmetric effective wake calculation performed in each 
pie segment. This assumption is least appropriate in the 0 
degree region of the wake. A more general effective wake 
calculation scheme may be required to obtain better results 
in this region. 

Results of M1T-PUF-3 cavitation performance of the 
AO-177 propeller are shown in Appendices III and IV. As an 
indicator of the influence of effective wake modifications 



to the 


n om i n a 1 


wa k e , 


MET-PUF-3 


wa s 


run 


twice, once with 


the 


n om i n a I 


wa k e as 


input 


and once 


with 


the 


effective wa k e 


(the 


mo re accurate 


bounda r y cond i t i 


on ) 


a s 


input. Appendix 


I I I 



shows cavitation performance given the nominal wake as 
input while Appendix IV shows cavitation performance given 
the effective wake as input. In general, the plots show 
that the cavitation extent is not changed much, although 
cavitation performance is better in the effective wake at 
the lower radii where, as has been shown, the effective 
wake modifications are the greatest. At the outer radix, 
where the effective wake modification is the least 
significant, the cavitation performance is least changed 
between the nominal wake and effective wake cases. 

Appendix V contains a graph comparing cavitation 
extent of MIT-PUF-3 output given nominal wake input against 
that given effective wake input. It plots the length of 



cavitation on the seventh radial cavitation panel (about 
the .8 propeller radius) non-dimensional ized on the local 
chord length as the blade rotates through 360 degrees. 
This graph indicates that the effective wake modification 
did not reduce cavitation extent on a blade, a result noted 
earlier. However, a larger contrast between the nominal 
wake and effective wake cases is shown in Appendix VI which 
graphs the cavitation volume ort a blade in the two cases. 
The volume is non-dimensional ized on the cube of the 
propeller radius. This graph shows a drop of approximately 
35% of peak cavitation volume. This change in cavity 
volume demonstrates the sensitivity of cavitation 
performance to small changes in the wake profile. This 
confirms that efforts to refine effective wake calculations 
were indeed justified. 

The effective wake cavitation performance indicated by 
MIT-PUF-3 compares reasonably well with experiments 
performed by the SSPA propeller tunnel. These experiments 
measured cavitation performance of the AO-177 propeller by 
photographing model tests. While no cavitation volume 
results were obtained, cavitation extent could be compared 
with MIT-PUF-3 oredict ions . Despite reductions of the 
effective wake calculations compared with nominal wake 
c a I c u I a t i o n s , MIT-PUF-3 still tended to predict greater 
cavitation than the experiments indicated, particularly at 
the middle radii. This seems to indicate that the 
effective wake modifications, while tending to the proper 
direction, were not large enough. This is consistent with 



the conclusions from comparing the two wakes that 



a mo r e 



general ef 
high wa k e 



fective wake scheme including rad 
regions is required. 



ial vorticity in 



CONCLUSIONS 



V . 



1 ) 

a k i s ymme t 
propel ter 

rt s t i 1 

modificat 
degrees . 
calculat i 
be r e qu i r 



Calculation of an effective wake using Huang's 
r ic method in pie-shaped segment with unsteady 
-induced velocities is a reasonable approximation. 
1 tends to underestimate the effective wake 
ion to the nominal wake in the vicinity of 0 
To correct this, a more general effective wake 
on scheme including radial vorticity effects may 
e d . 



2) The current effective wake calculation scheme can 
be used to predict cavitation pe-rf ormance for highly skewed 
propellers with 'up to seven blades. 



3) Cavitation 
small changes in 
improvements in 

signif icant ly mo r e 

predict ions . 



per f o rmanc e i s 

the wa k e prof i 
ef feet ive wa k e 
accurate cav 



1 



i 



indeed sensitive to 
e. Seemingly small 
calculations yield 
tation performance 
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V I . 



RECOMMENDATIONS 



1 ) 

including 



A mo re general 
the effects of 



effective walce calculation scheme, 
radial vorticity, is required. 



2) The computational time (and hence cost) of the 
unsteady field point velocity program (UFPV) can probably 
be cut without significant sacrifice of accuracy by 
decreasing the number of points calculated per radius. 
Appendix I shows that the variation of induced velocity 
with rotation is slight and smooth. Calculation of fewer 
points and estimation of the remaining points using 
harmonic analysis would probably not reduce 



accuracy . 



** APPENDIX I ** 



Norn ina 1 



, Induced, and Effective Velocities 



-21 - 




ooam/ch^a 3avn ivixy 



►4 

O 

>000 



1—4 1— I 1— « 



e *3 

O C Cjh 
2: M ud 



22 



NOM, IND, EFF VELOCITIES, .40 RADIUS 




ClOA^A/CtfM 2XYA TVIXV 



O 

>000 

«. / « • • 

> <fM * i ~* 



s 

o 

ss 



^ 
C ^ 
M W 



23 



NOM, IND, EFF VELOCITIES, .60 RADIUS 




C10A)A/C*DA 1VIXV 



pA 

o 

>000 



>* <t-< *<rH T— t 



E 

O 

2: 



-d 

C Cm 
M Cd 



24 



NOM, IND, EFF VELOCITIES, .70 RADIUS 




CIOAWCSM 3MVtt 1VIXV 



m3 

o 

>000 

"> — / • • • 

iM 



s 

o 



'S 



Cm 

Cm 

M 



25 



NOM, IND, EFF VELOCITIES, .80 RADIOS 




ClOAWCaM 2XVrt 1VIXV 



o 

>000 

V • • • 

» 1— t 



£ Tj 'h 
O c Cm 
S H W 



- 26 - 



NOM, IND, EFF VELOCITIES. .96 RADIUS. 




CIOAM/C^DA 1VIXV 



o 

>000 

s*/ • • « 

^ ^r-f 1—1 



o c 3 h 

X M 2 



27 



NOM, IND. EFF VELOCITIES, .98 RADIUS 



** APPENDIX II ** 

Nominal, Induced, and Effective Velocities, 
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** APPENDIX V ** 



Cavity Length, Nominal and Effective Valces 
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Cavity Vol ume , 



Nominal and Effective Wakes 
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56 



BLADE ANGLE 

NOMINAL VS. EFFECTIVE WAKE 



SSPA 



** A P PEND I XVII ** 



Experimental Propeller Cavitation Extent Plots 
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Nominal VaVre Diagram 
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Actual Velocities (i.e., V(VOL)-l.O) 
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Actual Velocities (l.e., V(V0L)=1.0) 
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MAIN PROGRAM IN FPUV. A PROGRAM TO CALCULATE 
UNSTEAOY VELOCITY HARMONICS AT FIELO POINTS 
SPECIFIEO INTERACTIVELY IN ROTERM. 



i ****** i 



UFPOOOIO 
UFP00O20 
UFP00030 
UFPOOO40 
UFP00050 
UFPOOO 6 O 
UFPOO070 
UFPOOO 8 O 
UFPOO09O 
UFPOO lOO 
UF PCX) 1 10 
UFPOO 1 20 
UF POO 1 30 
UFPOO 1 40 
UFPOO 1 50 
UFPOO 1 60 
UFPOO 1 70 
UFPOO 180 
UFPOO 1 90 
UFP00200 
UFP00210 
UFP00220 
UFP00230 
UFP00240 
UFP00250 
UFP00260 
UFP00270 
UFP00280 
UFP00290 
UFP00300 
UFP003 10 
UFP00320 
UFP00330 

ISTOP IS A FLAG SET INTERACTIVELY IN ROTERM TO INOICATE ENO OF PR0UFP0O34O 

UFP00350 

1ST 0P = O UFP00360 

UFPCO370 

ROFILE WILL RE AO OISX OATA FILE CREATEO BY PUF-2. NOTE THAT NSR IUFP0O380 
PUF-2 INPUT MUST BE DIVISIBLE BY THE NUMBER OF BLAOES. N8LA0E . UFP00390 

UFP00400 
UFP004 10 
UFP00420 
UFP00430 
UFP00440 
UFP00450 
UFP004 60 
UFP00470 
UFP00480 
UFP00490 
UFP00500 
UFP005 1 O 
UFP00520 
UFP00530 



*** COMMON TO CONTAIN GEOMETRIC PARAMETERS 

COMMON / GEOCOM / MM, MN. NN. NW . NTORL , X( 1 2 . 10) . Y ( 12 , 10) , Z ( 12 , 10) . 
& XTIP( lO, 1 1 ) , YTIP( 10, 1 1 ) , ZTIP( 10. 1 1 ) , XH( 101 . 7 ) , YH( 10 1 . 7 ) , 

& ZH( lOI .7 ) ,XW( 21 , 10) . YW( 2 1 . 10) . ZW(2 1 , 10) . RH 

*** COMMON TO CONTAIN SINGULARITY STRENGTHS 

COMMON / SINCOM / GT ( 10 , 10, 60 ) , G8 ( 10 , 9 . 60) . GTV( 10) , SB( 10 . 9 ) . GTV 
& . , GAM(9.60) , GMEAN(9) , 0KFW(2O) , HBWX , HBWY , HBWZ 

*** COMMON TO CONTAIN FIELO POINT DATA, NBLAOE , IOENT, NSR, ANO VSR 

COMMON / FPCOM / N8LA0E , NSR , XF P , RF P , TZF P , YB AR , Z8 AR . 

& OBLAOE , OTFP .-IOENT ( 1 8 ) . TFP ( 60 ) . VSR. I PLOT ,01 AM 

*** COMMON TO CONTAIN VELOCITY OATA 

COMMON / VELCOM / VIX ( 60 ) , VI T ( 60 ) . V I R ( 60 ) , U ( 60 , 3 ) , NSAMP 

COMMON TO CONTAIN HARMONICS COEFFICIENTS 

COMMON / HARCOM / A ( 1 5 , 3 ) , 8 ( 1 5 . 3 ) , AMP ( 1 5 . 3 ) , PH( 15 , 3 ) , NH 



CALL ROFILE 

*** INTERACTIVELY RE AO PARAMETERS FOR TERMINAL SESSION 



905 FORMAT ( / ' ENTER NBLAOE. NTORL. IPLOT. 7 / 
& ' FOR HELP, ENTER ZERO FOR NBLAOE. 7 ) 

30 WRI TE ( 6 , 905 ) 

REAO ( 5 , * ) NBLAOE, NTORL. IPLOT 
IF (NBLAOE. LT.O) GOTO 20 
IF (N8LA0E . EO.O) GOTO 90 
IF (IPLOT. EO. 2) CALL P I EOAT ( I STOP ) 

IF (ISTOP. EO. 1 ) GOTO 20 
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*** ROTERM GETS OATA INTERACTIVELY FROM TERMINAL ANO COMPUTES 
NECESSARY FIELD POINT GEOMETRY. 

10 CALL RDTERM( ISTDP ) 

*** IF OPERATOR INDICATES ENO OF PROGRAM. ISTOP WILL BE SET TO 1 
IF (ISTOP, GT.O) GDTD 20 

RUN THROUGH FPSTEP FOR EACH "TIME STEP." 

DD 100 K= 1 , NSR 

FPSTEP CALCULATED VELOCITY INDUCED AT FIELD POINT BY ONE BLADE 
AT "TIME STEP" K 

100 CALL FPSTEP ( K ) 

*** SUMVEL WILL SUM FIELD POINT VELOCITIES INOUCED BY EACH BLADE. 
CALL SUMVEL 



UFP00540 
UFP0O550 
UFP00560 
UFP00570 
UFP00530 
UFP0O590 
UFP00600 
UFP006 10 
UFP006 20 
UFP0O630 
UFP0O540 
UFPOO650 
UFPOO66O 
UFP00670 
UFP0O630 
UFPOO690 
UFP00700 
UFP007 10 
UFP00720 
UFP0O730 
UFPOO740 



UFP0O750 

FDUCAL WILL GENERATE HARMONIC COEFFICIENTS FOR INOUCED FIELD PO I NTUF P007 60 
VELOCITY AS PROPELLER ROTATES. * UFP0O770 

UFP0O730 

CALL FDUCAL UFP00790 

UFP00800 

**• PLDTVL PLOTS INOUCED FIELD POINT VELOCITIES AS A FUNCTION DF BL A0EUFP0O8 10 

ROTATION. UFP0O820 

UFP00830 

901 FORMAT ( ' AXIAL CODRO I NATE ' . 4X . F5 . 3/ ' RADIAL CODRD I NATE ' . 2X . F6 . 3/ UFPO0S4O 

& ' ANGULAR COORO I NATE ' , 2X , F5 . 1 /// ' AXIAL HARMON I C ' , 6X , F 8 . 5/ UFP00850 

& ' RADIAL HARMONIC' ,5X,F8.5/' TANGENTIAL HARMONIC '.F8.5) UFP0O860 

WRITE (6 ,901 ) XFP , RFP , TZFP , ( AMP( 1 ,U) , J= 1 . 3) UFP00870 

IF ( I PLOT . EQ . 1 ) CALL PLDTVL ( N8LADE , NSR . NH , TFP , TZFP , U , A . B , AMP , PH , UFP0O880 
& XFP, RFP, IDENT) UFP00890 

GDTD 10 UFP00900 

906 FORMAT ( / ' NBLADE < 0 => STOP' UFP00910 

& / 8X , '« 0 => HE LP ' / 8X , 9 > 0 => NUMBER OF BLADES'/ UFP0O920 

& ' NTORL = O TOTAL EFFECT ' /8X , ' * 1 *> THICKNESS ONLY'/ UFP00930 

& 8X , ' 3 2 => LOADING ONLY'/' I PLOT = 0 => DATA TD TERMINAL'/ UFP00940 

& 8X,'» 1 => PLOT -- NOT IMPLEMENTED 05/02/82 -- RL J ' UFP0O950 

& / 8X , ' = 2 *> CREATE PIEWAKE DATA FILE'/) UFP00960 

90 WR I TE ( 6 , 906 ) UFP0O970 

GOTO 30 UFP00980 

20 STOP UFP00990 

END UFPOIOOO 
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SUBROUTINE FPVEL(K,XF.YF,2F.VIX,VIY,VIZ) 

*** COMMON TO CONTAIN GEOMETRIC PARAMETERS 

COMMON / GEOCOM / MM , MN , NN, NW , NTORL . X ( 1 2 . 10 ) , Y ( 1 2 . 10) , 2 ( 1 2 , 10) , 
& XTIP( lO. 1 1 ) . YTIP( 10, 1 1 ) ,ZTIP( 10, 1 1 ) ,XH( 101 ,7) , YH( 101 ,7) , 

5 2H( lOI , 7 ) ,XW(2 1 . lO) . YW( 21 . 10) , 2W(21 , 10) ,RH 

*** COMMON TO CONTAIN SINGULARITY STRENGTHS ' • 

COMMON / SINCOM / GT ( 10 , 10. 60) . G8 ( 1 0. 9 . 60 ) . GTW( 1 0 ) . S3 ( 1 0 , 9 ) . GTV 

6 , GAM(9,60) , GME AN ( 9 ) t 0KFV(2O) , H8WX , HBWY , H8VZ 

*** COMMON TO CONTAIN FIELO POINT DATA, NBLAOE , IOENT , NSR , ANO VSR 

COMMON / FPCOM / NBLAOE , NSR , XFP , RFP , TZFP . YEAR , ZBAR t 
& OBLAOE , OTFP , IOENT ( IS ) . TFP ( 60 ) , VSR . I PLOT ,01 AM 

DIMENSION UT X ( 20 , 2 ) ,UTY(20,2) ,UTZ(20,2) ,U8X(20) .U8Y(20) ,UBZ(20) 
& , V0FWX(20,9) . V0FWY(2O, 9) ,V0FVZ(20,9) 

*** INITIALIZE VARIABLES 

MI0=MN/2 
NVFW=NW 
vJMPFW= 1 
WAKX -0 . 

WAKY=0 . 

WAKZ =0 . 

VIX=0.0 
V I Y = 0 . O 
VIZ=0.0 
DO 1 M*1,MM 
00 1 N-1.NN 

C CHOROWISE VORTICES ON 8LA0E 

I F ( NTORL . EO . 1 ) GO TO 8 

CALL VORSEG(XF, YF .ZF ,X(N,M) , Y(N, M) , Z(N.M) ,X(N+1 , M) , Y ( N+ 1 ,M) , 

1 Z(N+1 ,M) .CVX.CVY , CVZ. SX. SY , SZ.O) 

C SPANWISE VORTICES ANO SOURCES ON BLAOE 

8 CALL VORS E G ( XF ,YF,ZF,X(N,M),Y(N,M),Z(N t M) t X(N,M+l),Y(N f M-M), 

1 Z(N, M+ 1 ) , SVX , SVY , SVZ f SX , S Y , SZ , 1 ) 

S8NM- SB ( N , M ) 

G3NM=G8 ( N , M , K ) 

GTNM=GT ( N , M , K ) 

I F ( NTORL . EO . 0 ) GO TO 9 
IF(NTORL. EO. 1 ) G8NM=0.0 
IF(NTORL. EO. 1 ) GTNM=0.0 
I F (NTORL . EO . 2 ) SBNM = 0.0 

9 VIX~ VIX + G8NM* SVX+SBNM* SX + GTNM*CVX 
VI Y * V I Y + GBNM* SVY + SBNM* S Y + GTNM *CVY 
VIZ*VI Z + GBNM * SVZ+SBNM* SZ+GTNM*CVZ 

1 CONTINUE 

C SEPARATEO TIP VORTEX ON 8LA0E-- - - 



UF PO 1010 
UFP01020 
UF PO 1 030 
UF PO 1 040 
UFP01050 
UFP01060 
UFP01070 
UF PO 1 080 
UFP01090 
UFPO 1 lOO 
UFP01 1 10 
UFPO 1 120 
UFPO 1 130 
UFPO 1 140 
UFPO 1 150 
UFPO 1 160 
UFPO 1 170 
UFPO 1 180 
UFPO 1 190 
UFPO 1 200 
UFP012 10 
UFPO 1 220 
UFPO 1 230 
UFP01 240 
UFPO 1 250 
UFPO 1 260 
UFPO 1270 
UFPO 1 2 80 
UFPO 1 290 
UFPO 1 300 
UFPO 1310 
UFPO 1 320 
UFPO 1 330 
UFPO 1 340 
UFPO 1 350 
UFPO 1 360 
UFPO 1 370 
UFPO 1 380 
UFPO 1 390 
UFPO 1 400 
UFP014 10 
UFPO 1 4 20 
UFPO 1 4 30 
UFP01 440 
UFPO 1 4 50 
UFPO 1 4 60 
UF P0 1470 
UFPO 1 480 
UFPO 1 490 
UFP01 500 
UFPO 1510 
UFPO 1 520 
UFPO 1 530 



6 7 



non 



C 



I F ( NTORL . EO . 1 ) GO TO 10 

00 2 N= 1 ,NN 

CVX=0.0 

CVY =0 . O 

CVZ=0.0 

NTEM*NN-N+ 1 

00 3 LL- 1 , NTEM 

L 1 =LL+ 1 

CALL VORSEG(XF . YF , ZF , XT IP ( N , LL ) , YT I P ( N , LL ) , ZT I P ( N . LL ) 

* . XTI P ( N , L 1 ),YTIP(N,L1 ),ZTIP(N,L1 ),VX,VY,VZ,SX f SY.SZ.O) 

CVX=CVX+VX 
CVY *CVY+ VY 
3 CVZ=CVZ+VZ 

VIX=VIX+G3 ( N , MM , K ) *CVX 
VI Y = VI Y + GB ( N , MM t K ) *CV Y 
2 VI Z= VI Z+GB ( N , MM , K ) *CVZ 

TRANSITION WAKE 

WAKX -O . 

WAKY =0 . 

WAKZ=0 . 

00 35 N=1 , 100 

CALL VORSEG( XF , YF . ZF , XH( N , 1 ) , YH(N, 1 ) . ZH(N , 1 ) , XH( N+ 1 , 1 ) , YH( N+ 1 , 1 ) . 
1 ZH(N+1. 1) t VX,VY.VZ.SX,SY,SZ,0) 

WAKX = WAKX + VX 
WAKY =WAKY+VY 
35 WAKZ=WAKZ+VZ 

CALL VORSEG(XF . YF.ZF ,XH( 1 , 1 ) .0. .0. ,XH( 10 1 . 1 ) . O . , O . , 

& HBWX . HBWY . HBWZ ,SX,SY,SZ t O) 

HBWX= -HBWX 
HBWY = -HBWY 
HBWZ= -HBWZ 
11 00 71 M= 1 , MM 

I F ( M . EO . 1 ) GO TO 72 

00 73 N» 1 # NW 

UTX ( N , 1 )=UTX(N. 2) 

UTY ( N . 1 )«UTY(N,2) 

73 UTZ ( N , 1 )=UTZ(N.2) 

72 00 74 N* 1 ,NW 

CALL VORSEG( XF , YF . ZF , XW ( N , M ) . YW( N . M ) . ZW( N , M ) , XW ( N , M+ 1 ) , YW (N , M+ 1 ) . 
1 ZW(N t M+ 1 ) , USX ( N ) ,UBY(N) t UBZ(N) , SX.SY. SZ.O) 

IF(M.GT..1 ) GO TO 75 

CALL VORSEG (XF , YF , ZF t XW( N , M ) , YW(N.M) , ZW(N.M) ,XW(N+ 1 ,M) , YW(N+ 1 ,M) . 
1 ZW(N+1 ,M) , UTX ( N , 1) # UTY(N, l).UTZ(N, 1 ) , SX . SY , SZ , O ) 

75 CALL VORS EG( XF ,YF, ZF , XV( N t M+ 1 ) . YW( N , M+ 1 ) , ZW ( N , M+ 1 ) . XW ( N+ 1 , M+ 1 ) , 

1 YW ( N+ 1 , M+ 1 ) . ZW ( N+ 1 , M+ 1 ) , UTX ( N , 2 ) , UTY (N,2).UTZ(N,2),SX,SY.SZ,0) 

74 CONTINUE 

00 76 N» 1 ,NW 

*** COMPUTE INDUCED VELOCITY FROM ONE FORCESHOE AT THIS RAOIUS 

BUG=UTX(N,2) -UTX ( N , 1 )+UBX(N) 

CAT *UTY (N t 2)-UTY(N. 1 )+UBY(N) 

OOG=UTZ (N,2)-UTZ(N, 1 )+UBZ(N) 



UFP01540 
UFP01550 
UFP01560 
UF PO 1570 
UFP01580 
U F PO 1 590 
UF P0 1 600 
UF PO 16 10 
UFP01620 
UF P0 1630 
UFP01640 
UF P0 1 650 
UFPO 1660 
UFP01670 
UFPO 1 680 
UFPO 1 690 
•UFPO 1 700 
UFP017 10 
UF PO 1 7 20 
UFPO 1 730 
UFPO 1 7 40 
UFPO 1 750 
UFPO 1 7 60 
U F PO 17 70 
UFPO 1780 
UFPO 1 790 
UFPO 1 800 
UFP018 10 
U F PO 1 820 
UFPO 1 830 
UFPO 1840 
UF PO 1850 
UFPO 1 860 
UFPO 1 870 
UFPO 1880 
UFPO 1 890 
UFPO 1900 
UFP019 10 
UFP01920 
UFPO 1930 
UFPO 1 940 
UFPO 1950 
UFPO 1 960 
UFPO 1 970 
UFPO 1 980 
UFPO 1 990 
UFP02000 
UFP02010 
UFP02020 
UFP02030 
UFP02040 
UFP02050 
UFP02060 



63 



o o o o 



77 



73 



76 

92 



93 



71 

10 



IF(N.NE.NW) GO TO 77 UFP02070 

IF(M.NE.MID) GO TO 78 UFP02080 

UFP02090 

ALL TRANSITION WAKE VORTICES DISAPPEAR EXCEPT MID. WHICH TURNS UFP02100 
INTO ULTIMATE WAKE VORTICES (TIP ANO HUB) BY ASSUMING MEAN LOAD I NGUF P02 1 10 

UFP02 120 

BUG =BUG+WAKX + H8WX UFP02 130 

CAT =CAT +WAKY+HBWY UFP02 140 

DOG=DOG+WAKZ+HBWZ UFP02150 

GO TO 78 UFP02160 

BUG=BUG -UBX ( N+ 1 ) UFP02170 

CAT =CAT - UB Y ( N+ 1 ) UFP02 180 

DOG=DOG -UBZ ( N+ 1 ) UFP02190 

VOFWX (N,M)=BUG UFP02200 

VOFWY (N,M)=CAT UFP02210 

VOF WZ (N.M)=DOG UFP02220 

CONTINUE UFP02230 

DO 93 LN= 1 . NVFW UFP02240 

KKM=K-LN UFP02250 

IF (KKM.LE.O) KKM=KKM+NSR UFP02260 

VIX= VIX+VOFWX ( LN , M ) * ( GMEAN( M ) UFP02270 

1 +DKFW(LN)*(GAM(M.KKM)-GMEAN(M) ) ) UFP02280 

V I Y = VI Y + VOFWY ( LN , M ) * ( GMEAN( M ) UFP02290 

1 +OKFW( LN) * ( GAM( M , KKM ) -GME AN ( M ) ) ) UFP02300 

VI Z=VI Z+VOFWZ ( LN , M ) * ( GMEAN( M ) UFP023 10 

1 +OKFW( LN )*(GAM(M, KKM ) -GME AN ( M ) ) ) UFP02320 

CONTINUE ' UFP02330 

RETURN UFP02340 

END UFP02350 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C *** COMMON TD CDNTAIN GEOMETRIC PARAMETERS 
C 

CDMMON / GEOCDM / MM , MN . NN , NW . NTDRL , X ( 1 2 , 10 ) , Y ( 1 2 . 10 ) . Z ( 1 2 . 10 ) . 
& XTIP( 10, 1 1 ) . YTIP( 10, 1 1 ) , ZTIP( 10, 1 1 ) , XH( 101 .7) , YH( 101 , 7 ) . 

& ZH( 101 . 7) ,XW(21 . 10) . YV(21 . 10) ,ZW(21 , 10) ,RH 

CDMMON TD CONTAIN SINGULARITY STRENGTHS 

COMMDN / SINCOM / GT( 10 . 10, 60 ) ,GB ( 10 . 9 . 60 ) . GTW( 10) . SB ( 10 . 9 ) , GTV 
& , GAM (9, 60) , GME AN (9) . DKFW( 20 ) . HBWX f HBVY , HBWZ 

*** CDMMDN TO CONTAIN FIELD POINT DATA, N8LADE. IDENT. NSR t AND VSR 

CDMMON / FPCOM / N8LADE , NSR , XFP , RFP f TZFP . YBAR , ZBAR f 
& DBLADE . DTFP , IDENT ( 18) ,TFP(60) .VSR, I P LOT, D I AM 



SUBROUTINE RDFILE 



SUBROUTINE RDFILE READS DISK FILE14 DATA CREATED 
BY PUF-2 AND XAVES DATA NEEDED BY FPUV IN 
CDMMONS GEDCOM AND SINCDM. NOTE THAT NSR IN 
PUF-2 MUST BE DIVISIBLE BY THE NUMBER DF BLADES, 
NBLADE. 



NN => 

MM => 

NW => NUMBER OF VORTICES IN KEY WAKE 

NBLADE => NUMBER DF PRDP ELLER BLADES 

NVC *> NO DF CHORD VDRTICES, OTHER BLAOES 

NVS *> NO OF SPAN VDRTICES, OTHER BLADES 

NVW => NO OF WAKE VORTICES, OTHER WAKE 

NSR *> NUMBER DF TIME ST EP S/R EVOLUT I DN 

NWK => NUMBER DF RADII IN WAKE HARMONICS 



DIMENSION RZ( 10) , R ( 9 ) , CHDCP ( 9 ) ,CHORZ( 10) ,PHIB(90) 

DIMENSION XD(6,4,6),Y0(6,4,6),ZD(6,4,6) 

DIMENSION RWK ( 1 1 ) , AWA (16,11), BWA( 16.11), AWR ( 16,11), BWR (16,11), 
& AWT ( 16.11 ) ,BWT( 16,11) 

D I MENS I ON XWO (5.4,6), YWD (5.4,6), ZVO (5,4,6) 

DIMENSION SBO (5,3) ,G8C(5,3,60) , GAMWC (4,3, 60 ) ,GTC(5.4.60) . 

& GAMC(3, 60) , GMEANO ( 3 ) ,DK0W(4) 

READ (14) AU 
IF (AU.GE.2.0) STDP 

READ (14) I STDY , NBLADE , MM , NN , NW , NVC , NVS , NVW , NSR , NWK 
*** SET CONSTANT FDR GRID 



UFP02360 
UFP02370 
UFP02380 
UFP02390 
UFP0240O 
UFP024 10 
UFP02420 
UFP02430 
UFP02440 
UFP02450 
UFP02460 
UFP02470 
UFP02480 
UFP02490 
UFP02SOO 
UFP025 10 
UFP02520 
UFP02530 
UFP02540 
UFP02550 
UFP02560 
UFP02570 
UFP02580 
UFP02590 
UFP02600 
UFP02610 
UFP02620 
UFP02630 
UFP02640 
UFP02650 
UFP02660 
UFP02670 
UFP02680 
UFP02690 
UFP02700 
UFP027 10 
UFP02720 
UFP02730 
UFP02740 
UFP02750 
UFP02760 
UFP02770 
UFP02780 
UFP02790 
UFP02800 
UFP028 10 
UFP02820 
UFP02830 
UFP02840 
UFP02850 
UFP02860 
UFP02870 
UFP02830 



70 



*W(0M (DtnJkUKi 



NM=NN- 1 
MN- MM+ 1 
NP = NN + 1 
NZ=NW+1 
NVCP = NVC+ 1 
NVSP -NVS + 1 
NBLO = NBLADE - 1 
C 

C *** SET CONSTANTS FOR UNSTEADY PROBLEM 
C 

IF (ISTDY.EO.1) NSR = 1 
NVtyP=NVW+1 
NTB*NN*MM 
C 

C *-* READ KEY BLADE GEOMETRY 
C 

READ ( 14 ) ( (X(N,M) , N = 1 ,NP),M*1 , MN ) , 

2 ( ( Y ( N , M ) t N=1 , NP ) , M- 1 , MN ) , 

3 ( (Z(N,M) ,N*1 .NP) ,M«1 ,MN) , 

4 ‘ ( RZ ( M ) , M*1 , MN ) , 

5 (R(M) ,M*1 , MM) , 

6 (CHOCP(M) , M = 1 . MM) , 1 

7 (CHORZ(M) . M= 1 .MN) . 

8 (PHIB(K) , K- 1 , NTB ) 

C 

C *** READ BLADE GEOMETRY OF OTHER BLADES 
C 

IF ( NBLADE . LE. 1 ) GOTO 1 1 

RE AD ( 14) ( ( ( XD (N.M.K) , N = 1 . NVCP ) . M = 1 .NVSP ) .K = 1 , NBLD ) , 

2 ( ( ( YO( N.M.K) ,N=1 , NVCP ) . M= 1 .NVSP) ,K= 1 , NBLO ) , 

3 ( ( ( Z0( N.M.K) ,N= 1 .NVCP) . M* 1 .NVSP ) ,K=1 .NBLD) 

11 CONTINUE 

C 

C *** READ WAKE RADII AND WAKE HARMONICS 
C 

READ ( 14) (RWK(M) ,M= 1 , NWK ) 

READ ( 14) ( ( A WA ( I , M ) . 1*1 . 16) ,M=1 . NWK ) . 

( ( BWA ( I , M ) . I * 1 . 1 6 ) . M* 1 . NWK ) , 

( ( AWR (I.M).I-1,16),M=1, NWK ) . 

( ( BWR ( I. M). 1*1,16), M*1 .NWK) . 

( ( AWT ( I.M), 1*1,16), M*1 , NWK ) . 

( ( BWT ( I ,M) , 1=1 . 16) , M= 1 , NWK ) 

READ VARIOUS PROPELLER/WAKE CONSTANTS 

READ ( 14) WAKE , RPM . DI AM , RH , CDRAG , VSR , UR . ( I DENT ( N ) , N* 1 ,18) 

READ KEY WAKE GEOMETRY 

READ ( 14) ( ( XW ( N , M ) . N* 1 , NZ ) , M* 1 , MN ) , 

( ( YW( N ,M),N*1 ,NZ),M=1 ,MN), 

( (ZW(N,M) ,N=1 ,NZ) ,M*1 ,MN) . 

( (XH(N,K) ,N- 1 , 101 ) ,K*1 .NBLADE ) , 



UFP02890 
UFP02900 
UFP02910 
UFP02920 
UFP02930 
UFP02940 
UFP02950 
UFP02960 
UFP02970 
UFP02980 
UFP02990 
UFP03000 
UFP03010 
UFP03020 
UFP03030 
UFP03040 
UFP03050 
UFP03060 
UFP03070 
UFP03080 
UFP03090 
UFP03100 
UFP031 10 
UFP03120 
UFP03130 
UFP03 140 
UFP03 150 
•UFP03160 
UFP03 170 
UFP03180 
UFP03190 
UFP03200 
UFP03210 
UFP03220 
UFP03230 
UFP03240 
UFP03250 
UFP032S0 
UFP03270 
UFP03280 
UFP03290 
UFP03300 
UF P033 1 0 
UFP03320 
UFP03330 
UFP03340 
UFP03350 
UFP03360 
UFP03370 
UFP03380 
UFP03390 
UFP03400 
UFP034 10 
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c 

c 

c 



c 

c 

c 



c 

c 

c 



5 ( ( YH(N,K) ,N»1 , 101 ) ,K = 1 , N8LADE ) . 

6 ( (ZH(N.K) ,N-1 , 101 ) ,K* 1.NBLADE) . 

7 ( (XTIP(N, L) , L» 1 ,NP) ,N* 1 .NN) , 

8 ( (YTIP(N, L ) , L = 1 . NP ) ,N = 1 , NN) , 

9 ( ( ZTIP ( N , L ) , L* 1 , NP ) . N“ 1 , NN ) , 



*•* RE AO KEY 8LADE SINGULARITY STRENGTHS 

* ( (SB(N.M) , N- 1 , NN ) ,M= 1 t MM ) , 

1 ( ( (GB(N,M, KK ) ,N= 1 , NN ) , M = 1 , MM ) , KK= 1 , NS R ) , 

2 ( ( (GT(N.M.KK) ,N= 1 . NN ) . M = 1 , MN ) , KK- 1 ,NSR) . 

3 ( (GAM(M.KK) , M« 1 f MM ) . KK * 1 .NSP) 



*** OTHER WAKE GEOMETRY 



IF (NBLAOE.LE. 1 ) GOTO 12 

READ ( 14 ) ( ( (XWO( N.M . K) .N* 1 .NVWP ),M«1 .NVSP) ,K= 1 . NBLO ) , 

2 ( ( ( YWO (N,M,K) ,N= 1 ,NVWP ) ,M= 1 .NVSP ) ,K = 1 , NBLO ) , 

3 ( ( (2W0(N,M,K) ,N« 1 t NVWP) , M« 1 ,NVSP) .K= 1 .NBLO) . 

*** OTHER BLADE SINGULARITIES 



4 ( ( SBO (N.M), N= 1 . NVC ) ,M* 1 ,NVS) t 

5 ( ( (GBC(N.M.KK) , N= 1 , NVC ) , M= 1 , NVS ) ,KK» 1 ,NSR) . 

6 ( ( ( GAMWC (N.M, KK ) ,N= 1 . NVW ) , M= 1 ,NVS J.KK-1.NSR), 

7 ( ( (GTC(N,M,KK) ,N=1 , NVC ) , M= 1 , NVSP ) , KK- 1 ,NSR ) . 

8 ( (GAMC(M.KK) ,M=1 .NVS ) t KK=1 . NSR) , 

9 (GMEAN(M) f M=1.MM), 

* (GMEANO(M) , M= 1 , NVS ) , 

1 (DKFW(N) .N= 1 ,NW), 

2 (DKOW(N) , N= 1 .NVW) 

12 CONTINUE 

99 GOTO lOO 



C 

C &&& DIAGNOSTICS &&&&&&&&&&&&&&&&&&&&&&& , 

C 

901 FORMAT ( 1615) 

902 FORMAT ( 8F 10.5) 

903 FORMAT ( 1 8 A4 ) 

WR I T£ ( 1 . 902 ) AJ 

WRITE ( 1 , 901 ) ISTDY.NBLADE .MM.NN.NW.NVC.NVS .NVW.NSR. NWK 
WRITE ( 1 . 902) ((X(N,M),N-1 ,NP),M S 1 , MN ) . 

2 ( (Y(N.M) ,N«1 .NP) ,M*1 ,MN) , 

3 ( (2 (N.M) ,N« 1 .NP ) t M= 1 ,MN) , 

4 (RZ(M) t M*1 .MN) . 

5 ( R ( M ) , M= 1 , MM ) . 

6 (CHOCP(M) ,M= 1 .MM) , 

7 ( CHORZ ( M ) t M= 1 t MN ) f 

3 (PHIB(K) ,K* 1 , NTB ) 

IF (NBLADE.LE. 1 ) GOTO 21 

WRITE ( 1 .902) ( ( (XO ( N t M , K ) ,N-1 .NVCP) , M» 1 .NVSP) .K*1 .NBLO) . 

2 ( ( ( YO(N.M.K) ,N* 1 .NVCP ) . M s 1 .NVSP ) . K- 1 .NBLO ) . 

3 ( ( (ZO(N.M.K) ,N=1 .NVCP) .M= 1 ,NVSP) .K*1 .NBLO) 



UFP03420 
UFP03430 
UFP03440 
UFP03450 
UFP03460 
UFP03470 
UFP03480 
UFP03490 
UFP035OO 
UFP035 10 
UFP03520 
UFP03530 
UFP03540 
UFP03550 
UFP03560 
UFP03570 
UFP03580 
UFP03590 
UFP036OO 
UFP036 lO 
UFP03620 
UFP03630 
UFP03640 
UFP03650 
UFP03660 
UFP03670 
UFP03680 
UFP0369O 
UFP03700 
UFP037 10 
UFP03720 
UFP03730 
UFP03740 
UFP03750 
UFP03760 
UFP03770 
UFP03780 
UFP03790 
UFP03800 
UFP038 10 
UFP03820 
UFP03830 
UFP03840 
UFP03850 
UFP03860 
UFP03870 
UFP03880 
U FP03890 
UFP03900 
UFP03910 
UFP03920 
UFP03930 
UFP03940 



72 



21 CONTINUE 

WRITE (1,902) (RWK(M) ,M=1 , NWK ) 

WR ITE ( 1 . 902 ) ( ( AWA ( I , M ) , I = 1 , 16), M= 1 , NWK ) , 

3 ( ( BWA ( I , M ) , I = 1 , 1 6 ) , M= 1 , NWK ) , 

4 • ( ( AWR( I ,M) . I = 1 , 16) ,M= 1 .NWK) , 

5 ((BWR(I.M), 1=1, 16) ,M»1 .NWK). 

6 ( ( AWT( I ,M) . 1= 1 . 16 ) ,M= 1 .NWK) , 

7 ( (BWT( I ,M). 1= 1 . 16) ,M=1 .NWK) 

WR I TE ( 1 .902) WAKE.RPM.DIAM.RH.CDRAG.VSR.UR 
WRITE ( 1 , 903 ) (IDE NT (N) ,N=1 . 18) 

WRITE ( 1.902) ( ( XW( N ,M) ,N=1 ,NZ) ,M=1 , MN ) , 

2 ( (YW(N.M) ,N=1 ,NZ) ,M=1 ,MN) , 

3 ( ( ZW( N , M ) , N= 1 , NZ ) , M= 1 . MN ) , 

4 ( (XH(N.K) ,N= 1 . 10.1 ) ,K = 1 . NBLADE ) . 

5 ( ( YH( N , K ) , N= 1 , 101 ) ,K=1 .NBLADE) . 

.6 ( (ZH(N.K) .N= 1 . 101 ) ,K= 1 .NBLADE) , 

7 ( ( XT IP(N.L),L=1, NP ) . N= 1 , NN ) , 

8 ( (YTIP(N.L) ,L«1 ,NP) ,N« 1 .NN) , 

9 ( (ZTIP(N.L) ,L=1 .NP) ,N= 1 .NN) 

WRITE ( 1 , 902 ) ( ( SB(N.M) ,N= 1 ,NN) ,M=1 ,MM) , 

2 ( ( (GB(N.M.KK) ,N* 1 ,NN) ,M= 1 , MM ) ,KK= 1 ,NSR ) , 

3 ( ( (GT(N.M.KK) ,N= 1 ,NN) , M= 1 ,MN) ,KK=1 ,NSR ) , 

4 ( (GAM(M.KK) ,M=1 ,MM) ,KK= 1 .NSR) 

IF (NBLADE. LE. 1) GDTD 22 

WRITE ( 1 , 902) ( ( (XWO(N.M.K) ,N= 1 .NVWP) .M= 1 ,NVSP ) , K = 1 , NB LO ) , 

2 . ( ( ( YWD( N, M , K ) , N= 1 , NVWP ) . M= 1 , NVSP ) . K= 1 , NBLD ) , 

3 ( ( ( ZWO(N.M.K) , N= 1 .NVWP) ,M=1 .NVSP) , K= 1 . NBLO ) 

WR I TE ( 1 ,902) ( (SBD(N.M) ,N=1 .NVC) ,M=1 .NVS) 



2 

3 

4 

5 

6 

7 

8 
9 



22 

100 



CONTINUE 

RETURN 

END 



( ( (GBC(N.M.KK) ,N= 1 .NVC) ,M= 1 ,NVS ) ,KK= 1 .NSR) , 

( ( (GAMWC(N.M.KK) ,N= 1 .NVW) ,M= 1 .NVS ) ,KK= 1 ,NSR ) 
( ( (GTC(N.M.KK) ,N=1 , NVC ) .M= 1 .NVSP) ,KK« 1 .NSR ) . 
( (GAMC(M.KK) ,M=1 , NVS ) . KK= 1 . NSR ) . 

(GMEAN(M) ,M= 1 , MM ) . 

(GMEANO(M) ,M= 1 .NVS) . 

(DKFW(N) ,N=1 ,-NW) , 

(DKOW(N) ,N=1 .NVW) 



UFP03950 
UFP03960 
UFP03970 
UFP03980 
UFP03990 
UFP04000 
UFP04010 
UFP04020 
UFP04030 
UFP04040 
UFP04050 
UFP04060 
UFP04070 
UFP04080 
UFP04090 
UFP04 lOO 
UFP04 1 10 
UFP04 120 
UFP04 130 
UFP04 140 
UFP04 150 
UFP04 160 
UFP04 170 
UFP04 180 
UFP04 190 
UFP042OO 
UFP042 10 
UFPO4220 
UFP04230 
UFP04240 
UFP04250 
UFP04260 
UFP04270 
UFP04280 
UFP04290 
UFP04300 
UFP04310 
UFP04320 
UFP04330 
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SUBROUTINE ROTERM( I STOP ) 



UFP0434Q 



c 












UFP04350 


c 


* 








« 


UFP04360 


c 


* 


SUBROUTINE 


ROTERM RE AOS IN BLAOE ANO GEOMETRY 


UFP04370 


c 


« 


OATA 


INTERACTIVELY ANO COMPUTES NECESSARY FIELO - 


UFP04380 


c 


* 


POINT 


GEOMETRY. * 


UFP04390 


c 


* 








* 


UFP04400 


c 


* 


NBLAOE 


= > 


NUM8ER OF 8LA0ES * 


UFP044 10 


c 


* 


NTORL 


= 0 


= > 


CONS I OER THICKNESS ANO LOAOING * 


UFP04420 


c 


* 




* 1 


= > 


CONSIOER THICKNESS ONLY * 


UFP04430 


c 


* 




* 2 


a> 


CONS I OER LOAOING ONLY * 


UFP04440 


c 


* 


XFP 




= > 


AXIAL FIELO POINT COORDINATE, 


UFP04450 


c 


* 








POSITIVE OOWNSTRE AM * 


UFPO4460 


c 


* 


RFP 




*> 


RAO I AL FIELD POINT COORDINATE * 


UFP04470 


c 


* 


TZFP 




= > 


INI TAL FIELO POINT ANGLE IN OEG . 


UFP04430 


c 


* 


YBAR , 


ZBAR 




Y ANO Z FIELD POINT COORDINATES * 


UFP04490 


c 


* 


TH 




= > 


ANGLE FROM TIME STEP INPUT 


UFP04500 


c 


* 


YB, ZB 




= > 


Y , Z FIELO PT COORO, BEFORE TIME ST * 


UFP045 10 


c 

c 


* 








« 


UFP04520 

UFP04530 



C 

C 

c 



& 

& 



COMMON TO CONTAIN GEOMETRIC PARAMETERS 

COMMON / GEOCOM / MM . MN ? NN , NW , NTORL , X ( 1 2 . 1 0 ) . Y ( 1 2 , 10) . 2 ( 1 2 , 10 ) , 
XTIP< 10, 1 1 ) , YTIP( 10, 1 1 ) , 2TIP( 10. 1 1 ) ,XH( 101 .7 ) , YH( 101 . 7) . 

ZH( 101 , 7 ) , XW( 21 , 10) . YW( 2 1 , 10) , 2W( 21 , 10) , RH 

«*• COMMON TO CONTAIN FIELO POINT OATA, NBLAOE, IDENT , NSR , ANO VSR 

COMMON / FPCOM / NBLAOE , NSR , XFP , RFP , TZFP , YBAR f 2BAR , 

& OBLAOE , OTFP . IOENT( 18) ,TFP(60) , VSR . I PLOT. 01 AM 

901 FORMAT ( / ' ENTER XFP, RFP, TZFP . * / 

& ' FOR HELP, ENTER -11. FOR XFP.') 

10 WR I TE ( 6 , 90 1 ) 

READ ( 5 , *) XFP, RFP, TZFP 
IF (XFP. LT. -10. ) GOTO 90 
IF (TZFP. LT. 0.0) GOTO 70 
IF (RFP.GT.O. ) GOTO 60 
20 WR I TE ( 6 , 902 ) 

902 FORMAT (/ ' ENTER NBV , N . M . KK ' / ' FOR HELP, ENTER ZERO FOR NSW . ' ) 

RE AO ( 5 , * ) NBW , N , M , KK 

TH = 6 . 283185*FL0AT( KK - 1 )/FLOAT ( NSR ) 

IF (NBW.EO.O) GOTO 80 

*** FIELO POINT CENTEREO ON BLAOE GRIO. 

IF (NBW . EO. 1 ) CALL FPGRI 0 (X , Y , Z , N . M , KK , 1 2 ) 

*** FIELD POINT CENTERED ON TRANSITION WAKE GRIO. 

I F ( N8W .£0.2) CALL FPGR I 0 ( XW , YW , ZV , N . M , KK , 2 1 ) 



UFP04540 
UFP04550 
UFP04560 
UFP04570 
UFP04580 
UFP04590 
UFP04600 
UFP046 10 
UFP04620 
UFP04630 
UFP04640 
UFP04650 
UFP04660 
UFP04670 
UFP04680 
UFP04690 
UFP0470G 
UFP047 10 
UFP04720 
UFP04730 
UFP04740 
UFP04750 
UFP04760 
UFP04770 
UFP04780 
UFP04790 
UFP04800 
UFP048 10 
UFP04820 
UFP04830 
UFP04840 
UFP04850 
UFP04860 



- 7 ^ - 



60 

70 

80 

903 

90 

904 



1 => BLAOE WAKE GRID'/SX, 
N',7X,'=> CHOROWISE COOROINATE 



I F ( NBW . NE . 1 . AND . NBW . NE . 2 ) GOTO 10 
DTFP=360. 0/ FLOAT ( NSR ) 

RETURN 
I STOP 2 1 
RETURN 
WR I TE ( 6 . 903 ) 

FORMAT ( / ' NBW = O 2 > HELP'/5X,' 2 
& ' = 2 => TRANSITION WAKE GRIO'/' 

& ' M',7X,'=> SPANWISE COORDINATE'/' KK'*SX, 

& '=> TIME STEPS ROrATEO + 1'/) 

GOTO 20 

WR I TE ( 6 t 904 ) 

FORMAT ( / ' XFP <-lO=> HELP # /8X, 

' >- 10=> AXIAL FIELO POINT COOROINATE, POSITIVE OOWNSTREAM ' / 

< O => SELECT COOROINATE USING BLAOE OR WAKE GRIO'/ 

=> RADIAL FIELD POINT COORDINATE'/ 

< 0 => STOP ' / 8X , ' > O 2 > ANGULAR FILE POINT COOROINATE 



' RFP 
8X , ' > O 
' T2FP 
/) 

GOTO 10 
ENO 



UFP04870 
UFP04880 
UFP04890 
UFP04900 
UFP049 10 
UFP04920 
UFP04930 
/UFP04940 
UFP04950 
UFP04960 
UFP04970 
UFP04980 
UFP04990 
UFPOSOOO 
UFP050 10 
UFP05020 
UFP05030 
UFP05040 
UFP05050 
UFP05060 
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SUBROUTINE FPGRIO(XX ,YY.ZZ.N,M.KK.NDIM) 

C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C *** CDMMDN TO CONTAIN FIELD POINT DATA, NBLAOE , I DENT , NSR , ANO VSR 
C 

CDMMDN / FPCDM / NBLADE . NSR , XFP , RFP , TZFP . YBAR , ZBAR , 

& OB LAO E ,DTFP, IDE NT ( 18) ,TFP( 60) , VSR , I PLOT , DI AM . 

C 

01 MENS I ON XX(NDIM, 10) .YY(NDIM, 10) ,ZZ(NOIM, 10) 

*** COMPUTE X , Y , Z CDDRO INATES IN ZERO ANGLE POSITION 

TH=6 .283185*FLDAT( KK- 1 )/FLDAT(NSR) 

THDEG = TH* 360 . /6 . 283 185 
OBLAOE=360.0/FLDAT (NBLAOE) 

XFP=0. 25* ( XX ( N , M ) + XX(N+1,M) + XX(N,M+l) + XX(N+1,M+1)) 

Y8 =0. 25* ( YY ( N , M ) + YY(N+1,M) + YY(N,M+1) + YY(N+1,M+1)) 

ZB =0 . 25* ( ZZ ( N , M ) + ZZ(N+1.M) + ZZ(N,M+l) + ZZ(N+1,M+1)) 

*** CONVERT TD CYLINDRICAL COORDINATES 

RFP=SQRT(Y8*Y3+Z8*ZB) 

TZFP- AT AN2 ( ZB , YB )/l . 7453293E-02 
TZFP=TZFP+THDEG 

*** BRANCH IF GENERATING PIEWAKE OATA 
IF ( I PLOT . EO . 2 ) GOTO 200 

TZFP = 1 80. 0/FLDAT( NBLAOE ) - AT AN2 ( ZB . YB ) / 1 . 7453293E-02 

*** KEEP TZFP BETWEEN 0 AND DNE BLADE ANGLE 

OD 100 K= 1 .NBLAOE 
IF(TZFP.LT.O.O) TZFP = TZFP+OB LAO E 
I F ( TZFP . GT . DBLAOE ) TZFP^TZFP -DBLAOE 
lOO CONTINUE 
C 



SUBROUTINE GRID CALCULATES X,R, AND THETA FIELD 
PD I NT GEDMETRY COMPONENTS GIVEN GRIO VECTORS 
X , Y , ANO Z FOR BLAOE OR WAKE. GRIO CDDRO I NATE S 
N AND M, AND TIME STEP KK . 

XX.YY.ZZ => GRIO GEDMETRY VECTORS 

N.M *> GRIO INDICES SELECTED FOR VELOCITY 

KK => TIME STEP. NUMBER OF ANGLE STEPS + 1 

XFP, YBAR, ZBAR , TZFP, RFP 

=> X.Y.Z, THETA, RADIAL FIELD PT CODRO 
OBLADE => ANGLE BETWEEN BLAOES 

TH => ANGLE OF TIME STEP 

NO I M => FIRST 0 1 MENS I ON DF GRIDS PASSED 



UFP05070 
UFP05080 
UFP05090 
UFP05100 
UFP051 10 
UFP05 120 
UFP05130 
UFP05 140 
UFP05150 
UFP05160 
• UFP05 170 
UFP05 180 
UFP05190 
UFP05200 
UFP052 10 
UFP05220 
UFP05230 
UFP05240 
UFP05250 
UFP05260 
UFP05270 
UFP05280 
UFP05290 
• UFP05300 
UFPO5310 
UFP05320 
UFP05330 
UFP05340 
UFP05350 
UFP05360 
UFP05370 
UFP05380 
UFP05390 
UFP054Q0 
UFP05410 
UFP05420 
UFP05430 
UFP05440 
UFP05450 
UFP05460 
UFP05470 
UFP05480 
UFP05490 
UFP05500 
UFP05510 
UFP05520 
UFP05530 
UFP05540 
UFP05550 
UFP05560 
UFP05570 
UFP05580 
UFP05590 
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*«* ROTATE BLADE FOR PIEWAKE DATA UFP05600 

UFPC56 10 

200 YBAR=YB»COS(TH) - Z3*SIN(TH) UFP05620 

ZBAR=YB»SIN(TH) + ZB*COS(TH) UFP05630 

300 RETURN UFP05640 

END UFP05650 
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C 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C *** COMMON TO CONTAIN GEOMETRIC PARAMETERS 
C 

COMMON / GEOCOM / MM , MN, NN , NW. NTORL , X( 12 . 10) , Y ( 1 2 , 10) , Z ( 12 , 10) , 
& XTIP( lO. 1 1 ) . YTIP( 10. 1 1 ) ,ZTIP( 10. 1 1 ) ,XH( 101 . 7) . YH( 101 ,7 ) , 

& ZH( 101 ,7) ,XW(21 , 10) . YW(21 t 10) ,ZW(21 , 10),RH 

*** COMMON TO CONTAIN SINGULARITY STRENGTHS 

COMMON / SINCOM / GT ( 10, 10 , 60 ) . G8 ( 10,9,60). GTW( 10), SB (10, 9), GTV 
& , GAM ( 9 , 60 ) , GMEAN( 9 ) , 0KFW(20) , HBWX , HBWY , H8WZ 

COMMON TO CONTAIN FIELD POINT OAT A . NBLAOE , I CENT , NSR , ANO VSR 

COMMON / FPCOM / NBLAOE , NSR . XFP , RFP , TZFP , YBAR , ZBAR , 

& OBLAOE , OTFP , IOENT( 18) , TFP ( 60 ) , VSR , I PLOT. 01 AM 

*** COMMON TO CONTAIN VELOCITY OATA 

COMMON / VELCOM / V I X ( 60 ) . VI T ( 60 ) . VI R ( 60 ) . U( 60 , 3 ) . NSAMP 

*** CALCULATE CARTESIAN FIELO POINT COOROINATES TO FEEO FPVEL 

TFP(K)*FLOAT(K- 1 )*OTFP-TZFP 
TRAD“TFP(K)*1 . 7 4 53 29 3 E -02 
COSFP=COS(TRAO) 

SINFP s SIN(TRAO) 

YFP 3 RFP*COSFP 
ZFP = RFP* SINFP 

FPVEL WILL CALCULATE VIX, VIY. ANO VIZ FOR FIELO POINT 

CALL FPVEL(K,XFP. YFP. ZFP. VIX(K) t VIY. VIZ) 

VIR(K)=VIY * COS FP + VI Z*SINFP 
VIT ( K ) = VIZ *COSFP - VIY * SINFP 
RETURN 
END 



SUBROUTINE FPSTEP(K) 



FPSTEP IS A SUBROUTINE TO CALCULATE THE FIELO 
POINT VELOCITIES INOUCEO BY ONE BLAOE AT “TIME 
STEP*. K. THE "OIRTY WORK “ IS OONE IN 
SUBROUTINE FPVEL. 

TFP *> ANGULAR FIELD POINT COORDINATE. OEG . 

TR AO => ANGULAR FIELO POINT COOROINATE. RAO. 

YFP. ZFP => Y ANO Z FIELO POINT COOROINATES 

VIX. VIY. VIZ.VIR. VIT => INOUCEO VELOCITY COMPONENT 



UFP05660 
UFP05670 
UFP05680 
UFP05690 
UFP05700 
UFP057 10 
UFP05720 
UFP05730 
UFP05740 
UFP05750 
UFP05760 
UFP05770 
UFP05780 
UFP05790 
UFP05800 
UFP058 10 
UFP05820 
UFP05830 
UFP05840 
UFP05850 
UFP05860 
UFP05870 
UFP05S80 
UFP05890 
UFP05900 
UFP05910 
UFP05920 
UFP05930 
UFP05940 
UFP05950 
UFP05960 
UFP05970 
UFP05980 
UFP05990 
UFP06000 
UFP060 10 
UFP06020 
UFP06030 
UFP06040 
UFP06050 
UFP06060 
UFP06070 
UFP06080 
UFP06090 
UFP06 100 
UFP06 1 10 
UFP06 120 
UFP06 130 
UFP06 140 
UFP06 1 50 
UFP06160 
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C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C *•* COMMON TO CONTAIN FIELO POINT DATA, NBLAOE , IDENT, NSR . ANO VSR 
C 

COMMON / FPCOM / NBLAOE t NSR ,’XFP . RFP , TZFP , YBAR , ZBAR , 

& OBLAOE,OTFP, IDENT ( 13),TFP(60) , VSR . I PLOT , 01 AM 

*** COMMON TO CONTAIN VELOCITY OATA 

COMMON / VELCOM / VIX ( 60 ) . VI T ( SO) , VI R ( 60 ) . U( €0 , 3 ) , NSAMP 

*** AT ANY GIVEN REAL TIME, EACH BLAOE CAN BE CONSIDERED TO HAVE 
TRAVELLED THROUGH A CERTAIN NUMBER OF "TIME STEPS * . 

* BY SUMMING VIX, VIR, ANO VIT AT BLAOE RATE, WE CAN INCLUDE 
THE VELOCITIES INDUCED BY EACH BLAOE INTO ONE VELOCITY 
INOUCED BY THE ENTIRE PROPELLER. 

NS AMP = NSR /NBLAOE 
DO 20 N= 1 , NSAMP 
00 10 .3 

10 U ( N , J ) =0 . O 

DO 20 K* 1 .NBLAOE 
L = N+(K-1 ) *NS AMP 
U(N, 1 )=U(N, 1 )+VIX(L) 

U(N,2)®U(N.2)+VIR(L) 

20 U(N,3)=U(N,3)+VIT(L) 

*** NOW CONVERT FROM REFERENCE VELOCITY TO SHIP VELOCITY 

00 30 N= 1 .NSAMP 
00 30 J=1 ,3 
30 U(N,U)®U(N,U)/VSR 
RETURN 
ENO 



SUBROUTINE SUMVEL 



SUBROUTINE SUMVEL SUMS VELOCITIES FROM EACH 
BLAOE, CALCULATED IN FPSTEP, TO PRODUCE THE 
TOTAL PROPELLER-INDUCED FIELO POINT VELOCITY. 



U(N, J) 



NSAMP 



VIT 


= > 




= > 


J-0 


= > 


J® 1 


= > 


J® 2 


= > 




= > 



VELOCITY COMPONENTS INOUCED BY 
ONE BLAOE AT "TIME STEP" K 
VELOCITY INOUCED 8Y ALL BLAOES 
AT "TIME STEP" N 
AXIAL COMPONENT 
RAOIAL COMPONENT 
TANGENTIAL COMPONENT 
NUMBER OF "TIME STEPS" IN ONE 
BLAOE ANGLE (= NUMBER OF 
U VELOCITIES COMPUTED.) 



UFP06 170 
UFP06 180 
UFP06 190 
UFP06200 
UFP062 10 
UFP06220 
UFP06230 
UFP06240 
UFP06250 
UFP06260 
UFP06270 
UFP06280 
UFP06290 
UFP06300 
UFP063 10 
UFP06320 
UFP06330 
UFP06340 
UFP06350 
UFP06360 
UFP06370 
UFP06380 
UFP06390 
UFP06400 
UFP064 10 
UFP06420 
UFP06430 
UFPQ6440 
UFP06450 
UFP06460 
UFP06470 
UFP06480 
UFP06490 
UFP06500 
UFP065 10 
UFP06520 
UFP06530 
UFP06540 
UFP06550 
UFP06560 
UFP06570 
UFP06580 
UFP06590 
UFP06600 
UFP066 10 
UFP06620 
UFP06630 
UFP06640 
UFP06650 
UFP06660 
UFP06670 
UFP06680 
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C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C *** COMMON TO CONTAIN FIELO POINT OATA , NBLADE . IDENT, NSR , ANO VSR 
C 

COMMON / FP COM / NBLAOE .NSR.XFP.RFP, T2FP, YBAR, 2BAR, 

& OBLAOE .OTFP, IOENT( 13) ,TFP(60) .VSR. I PLOT, 01 AM 

*** COMMON TO CONTAIN VELOCITY DATA 

COMMON / VELCOM / VIX ( 60) . VIT ( 60 ) . VI R ( 60) . U ( 60 . 3 ) . NS AMP 

*** COMMON TO CONTAIN HARMONICS COEFFICIENTS 

COMMON / HARCOM / A ( 1 5 , 3 ) , B ( 1 5 , 3 ) , AMP ( 1 5 , 3 ) , PH( 1 5 . 3 ) . NH 

*** FIRST CALCULATE A ANO B 

NFB=NSR/NBLAOE 
NHMAX = NF B/ 2+1 
NH=MINO( NHMAX , 15) 

STEP=6.283185/NFB 
ANFB=2 . O/NFB 
BNFB * 1 .O/NFB 
00 11 K s 1 , NH 
DO 10 J*1,3 
A ( K , J ) =0. O 
10 B ( K , J ) -O . O 

DO 9 N*1 , NFB 
T * ( (N-1 ) * ( K - 1 ) ) * ST E P 
ST*SIN( T ) 



SUBROUTINE FOUCAL 



SUBROUTINE FOUCAL TAKES THE PROPELLER INOUCED 
VELOCITIES CALCULATEO IN SUMVEL AND PERFORMS 
A HARMONIC ANALYSIS ON THEM. IT CALCULATES 
FOURIER COEFFICIENTS FOR BOTH COSINE (A) ANO 
SINE (B) ANO ALSO COSINE AMPLITUOE (AMP) AND 
PHASE ANGLE (PH). 



NH 

A(N,J) ,B(N,J) 



AMP ( N , J ) 

PH(N, J) 
NFB 



=> NUMBER OF HARMONICS CALCULATED 
*> FOURIER COEFFICIENTS 

A => COSINE COEFFICIENTS 
B =* SINE COEFFICIENTS 
N -> HARMONIC NUMBER UP TO 15 
J - 1 => AXIAL VELOCITY 

=» 2 *> RAO I A L VELOCITY 
* 3 => TANGENTIAL VELOCITY 
=> FOURIER COEFFICIENT AMPLITUDE 
N AND J SUBSCRIPTS AS FOR A, B 
«> PHASE ANGLE ASSOCIATED WITH AMP 
=> NUMBER OF FIELO POINTS / BLAOE 



UFP06690 
UFP06700 
UFP06710 
UFP06720 
UFP06730 
UFP06740 
UFP06750 
UFP06760 
UFP06770 
UFP06780 
UFP06790 
UFP0680O 
UFPQ68 10 
UFP06820 
UFP06830 
UFP06840 
UFP06850 
UFP06860 
UFP06870 
UFPO688O 
UFP06S90 
UFP06900 
UFP069 10 
UFP06920 
UFP06930 
UFP06940 
UFP06950 
UFP06960 
UFP06970 
UFP06980 
UFP06990 
UFP07000 
UFP070 10 
UFP07020 
UFP07030 
UFP07040 
UFP07050 
UFP07060 
UFP07070 
UFP07080 
UFP07090 
UFP07 100 
UFP071 10 
UFP07120 
UFP07130 
UFP07 140 
UFP07 150 
UFP07160 
UFP07170 
UFP07 180 
UFP07 190 
UFP07200 
UFP072 10 
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non 



CT=COS(T) 

00 9 J s 1 , 3 

A(K,d) s A(K,d)+U(N f d) *CT 
9 B(K t J)=B(K,J)+U(N,J)*ST 

*** CONVERT A AND B INTO AMP AND PH 

DO 11 d= 1 , 3 

IF(K.EO.I) A( 1 . J)=A( 1 , J)-BNFB 
IF(K. EO.NHMAX) GO TO 12 
IF(K.GT.I) A(K f d) s A(K f d) * ANFB 
IF(K.GT.I) B(K,d) s B(K f d)*ANFB 
GO TO 13 

12 A(K f d)*A(K, d)*5NFB 
8 ( K , d ) *0 . O 

13 AMP (K,d)=SQRT(A(K f d)**2+B(K,d)**2) 

PH(K, d)=ATAN2(A(K,d) , B(K,d) )*57.296 
IF(PH(K,d) -LT.O.O) PH( K f d) =PH(K f d) +360.0 
PH(K f d ) *PH(K, d )/NBLADE 

PH( 1,d)«0.0 
AMP ( 1 ,d)-A( 1 ,d) 

11 CONTINUE 
RETURN 
END 



UFP07220 
UFP07230 
UFP07240 
UFP07250 
UFP07260 
UFP07270 
UFP07280 
UFP07290 
UFP07300 
UFP073 10 
UFP07320 
UFP07330 
UFP0734O 
UFP07350 
UFP07360 
UFP07370 
UFP07380 
UFP07390 
UFP07400 
UFP074 lO 
UFP07420 
UFP07430 
UFP07440 
UFP07450 
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C 
C 
C 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C *** CDMMON TD CONTAIN GEOMETRIC PARAMETERS 
C 

COMMON / GEDCDM / MM , MN, NN # NW , NTORL , X ( 1 2 , 1 0 ) . Y ( 1 2 , 10 ) , Z ( 1 2 , 10 ) , 
& XTIP( 10 , 11),YTIP( 10,1l),ZTIP( 10, 11), XH( 101,7), YH( 101,7), 

& ZH( 101 f 7) ,XV(21 t 10) , YW(21 , 10) . ZW(21 . 10) , RH 

*** COMMON TD CONTAIN SINGULARITY STRENGTHS 

COMMON / SINCDM / GT( 10, 10. 60) , G8( 10, 9 . 60) , GTW( 10) . S3( 10 . 9) , GTV 
& , GAM (9,60). GMEAN( 9 ) , DKFW( 20 ) , H8WX , HBWY , HBWZ 

**- CDMMDN TD CONTAIN FIELD POINT DATA, NBLADE , IDENT, NSR , ANO VSR 

COMMON / FPCDM / NBLAOE , NSR , XFP . RF P , TZFP , YB AR , ZB AR , 

& DBLADE.DTFP, I DENT ( 18) ,TFP(60) , VSR, I PLOT . D I AM 

*** CDMMON TD CONTAIN VELOCITY DATA 

COMMON / VELCOM / VIX ( 60 ) , VI T ( 60 ) , VI R ( 60 ) . U( 60 , 3 ) , NS AMP 

*** CDMMON TD CONTAIN HARMONICS COEFFICIENTS 



SUBROUTINE P I EDAT ( I STOP ) 



SUBROUTINE PIEOAT GENERATES DATA FROM INPUT TO 
AXIAL VELOCITY DUTPUT TO BE USED AS INPUT FOR 
UNSTEADY PIEVAKE. IT USES BLADE GRID GEOMETRY 
TD CALCULATE NSR FIELD POINT AXIAL T IME -AVERAGED 
VELOCITIES AT FIVE RADII ALONG THE LEADING EDGE. 
IT THEN CALCULATES NSR FIELD POINT AXIAL TIME- 
AVERAGED VELOCITIES AT 1.1 TIMES PROPELLER RADIUS 

N => CHORDWISE GRID COORDINATE 

M => SPANVISE GRID COORDINATE 

K.KK => ANGLE INCREMENTS + 1 

I STDP* 1 => STOP FPUV UPON RETURN TD MAIN PROG 

NSR => NUMBER TIME STEPS IN REVOLUTION 

NR *> RADIUS NUMBER 

XFP, YBAR,ZBAR,RF P. TZFP 

=> X,Y,Z, RADIAL, THETA FIELD PT COORD 
THETA(1-5)=> ANGLES 5 TIME STEPS BEFORE ZERO 
BLADE ANGLE 

THETA ( 6 - NSR+5) => ANGLES IN ONE BLADE REV. 

THETA ( NSR+6 - NSR+20) => ANGLES 15 TIME STEPS 
AFTER ZERO BLAO ANGLE. 

VX ( K ) => AXIAL VELDCIT? AT ANGLE THETA(K) 

UA(KK,M) => SMOOTH AXIAL VELOCITY AT RADIUS M, 
ANGLE KK 



UFP07460 
UFP07470 
UFP07480 
UFP07490 
UFP075O0 
UFP075 10 
UFP07520 
UFP07530 
UFPQ7540 
UFP07550 
UFP07560 
UFP07570 
UFP07580 
UFP07590 
UFP07600 
UFP076 10 
UFP07620 
UFP07630 
UFP07640 
UFP07650 
UFP07660 
UFP07670 
UFP07680 
UFP07690 
UFP077OO 
UFP077 10 
UFP07720 
UFP07730 
UFP07740 
UFPO7750 
UFP07760 
UFP07770 
UFP07780 
UFP07790 
UFP07800 
UFP078 10 
UFP07820 
UFP07830 
UFP07840 
UFP07350 
UFP07860 
UFP07870 
UFP07880 
UFP07890 
UFPQ7 900 
UFP07910 
UFP07920 
UFP07930 
UFP0794O 
UFP07950 
UFP07960 
UFP07970 
UFP07980 
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c 

COMMON / HARCOM / A ( 1 5 . 3 ) . B ( 1 5 , 3 ) . AMP ( 1 5 . 3 ) . PH ( 1 5 . 3 ) . NH 
C 
C 

DIMENSION THETA (80) . VX ( 80 ) . UA ( 60 . 9 ) . RWK(9) 

N= 1 
NR* 1 
C 

C *** CALCULATE VELOCITIES AT FIVE BLAOE GRID RADII 
C 

102 00 300 M=3,9,2 

NR=NR+ 1 
C 

C *** CALCULATE VELOCITIES AT NSR ANGLES FOR EACH RADIUS 
C 

00 200 KK = 1 t NSR 
C 

C **• FIELD POINT CENTERED ON BLADE GRIO 
C 

CALL FPGR IO(X,Y.Z,N,M, KK ,12) 

IF (M.E0.3) X 1 = XF P 
KK 1 5 = KK+ 1 5 
THETA ( KK 15)=TZFP 
RWK ( NR ) =RF P 
DTFP=360. /FLOAT (NSR) 

C 

C *** CALCULATE AXIAL VELOCITY FOR FIELO POINT AT RAOIUS RWK ( NR ) ANO 
C *** ANGLE THET A ( KK 15). STORE THIS VELOCITY IN VX . 

C 

200 CALL AMPC AL ( VX ( 16) , KK ) 

NSR 1 5=NSR+ 1 5 
C 

C *** STORE SMOOTHED VELOCITIES IN NR -TH COLUMN OF UA 
C 

300 CALL SMOOTH( NSR. NR .VX.THETA.UA) 

C 

C *** CALCULATE NSR VELOCITIES AT PROPELLER HUB 
C 

NRHUB* 1 
XFP = 1 . 5*X 1 
RWK( NRHUB )=RH- 
C 

C *** STEP THROUGH NSR ANGLES AT HUB 
C 

DO lOO KK* 1 .NSR 
KK15=KK+15 

THETA (KK 15 ) =360. * FLOAT ( KK- 1 )/FLOAT(NSR) 

TZFP = THET A ( KK 1 5 ) 
lOO CALL AMPCAL(VX( 16) . KK ) 

NSR 1 5=NSR+ 1 5 
C 

C *** NOW STORE SMOOTHED HUB VELOCITIES IN FIRST COLUMN OF UA 
C 



UFP07990 
UFP08000 
UFP080 10 
UFP08020 
UFP08030 
UFP08040 
UFPC8050 
UFP08060 
UFP08070 
UFP08080 
UFP08090 
UFP08 100 
UFP08 1 10 
UFP08 120 
UFP08 130 
UFP08 140 
UFP08 150 
UFP08 160 
UFPC8 170 
UFP08180 
UFP08 190 
UFP082OO 
UFP08210 
UFP08220 
UFP08230 
UFP08240 
UFPO8250 
UFP08260 
UFP08270 
UFP082S0 
UFP08290 
UFP08300 
UFP083 10 
UFP08320 
UFP08330 
UFP08340 
UFP08350 
UFP08360 
UFP08370 
UFP08380 
UFP08390 
UFP08400 
UFP08410 
UFP08420 
UFP08430 
UFP08440 
UFP08450 
UFP08460 
UFP08470 
UFP08480 
UFP08490 
UFP08500 
UFP08510 



ooo nooooo non 



CALL SMOOTH ( NSR , NRHU8 , VX . THETA , UA ) 

C 

c NOW CALCULATE VELOCITIES AT A NEXT RADII OUTSIDE THE PROPELLER 

C *»» DISK. USE XFP OF PROPELLER GRID CASE OF N» 1 . M*9 . 

C 

C &&& DIAGNOSTICS &&&&&&&&&&&&&&&&& 

C 

C GOTO 101 

c 

C *** CALCULATE NSR VELOCITIES AT 1.1, 1.4, AND 1.7 PROP RADII 
C 

DO 500 KKK = 1 ,7,3 

RFP = 1 .0+FL0AT(KKK)/l0. 

NR=NR+ 1 
RWK ( NR ) =RFP 

*** CALCULATE NSR VELOCITIES AT THESE RADII 

DO 400 KK- 1 , NSR 
KK15=KK+15 

THETA (KK 15 ) = 360. * FLOAT ( KK - 1 )/ FLOAT (NSR) 

TZF P -THETA ( KK 1 5 ) 

400 CALL AMPCAL(VX( 16) ,KK) 

*** NOW STORE SMOOTHED VELOCITIES IN NR-TH COLUMN OF UA 
500 CALL SMOOTH(NSR. NR.VX.THETA.UA) 

*** WR I TE ( OUTPUT IN FORMAT FOR UNSTEADY PIEWAKE 

901 FORMAT ( 8 F 10 . 5 ) 

902 FORMAT ( 1615) 

101 WR I TE ( 3 , 902 ) NR 

WR I TE ( 8 , 90 1 ) ( RWK ( M ) , M 3 1 , NR ) 

WRITE (8 ,901 ) ' ( ( UA ( KK , M ) , KK= 1 .NSR) ,M=1 . NR ) 

*** STOP FPUV UPON RETURN TO MAIN PROGRAM 

I STOP 3 1 
RETURN x 

END 



UFP08520 
UFP08530 
UFP08540 
UFP08550 
UFP08560 
UFP08570 
UFP08580 
UFP08590 
UFP086O0 
UFP08610 
UFP08620 
UFP08630 
UFP08640 
UFP08650 
UFP08660 
UFP08680 
UFP08690 
UFP087OO 
UFP087 10 
UFP08720 
UFP08730 
UFP08740 
UFP08750 
UFP08760 
UFP08770 
UFP087 80 
UFP08790 
UFP08300 
UFP08810 
UFP08S20 
UFP08830 
UFP08840 
UFP08850 
UFP08860 
UFP08870 
UFP08380 
UFP08890 
UFP0890O 
UFP08910 
UFP08920 
UFP08930 
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SUBROUTINE AMPCA l ( VX , KK ) UFP08940 

C **.,,«**»,,«.,,«**»***.***.*««..»***«*.,.«*«»*».«.*,*, UFP08950 

C * * UFP08960 

C * SUBROUTINE AMPCAL WILL TAKE FIELO POINT GEOMETRY * UFP08970 

C * CALCULATED IN PIEOAT ANO CALCULATE THE INDUCEO - UFP08980 

C * TIME-AVERAGED AXIAL VELOCITY. IT WILL STORE THAT* UFP08990 

C * VELOCITY IN VX(KK). * UFP09000 

C * * UFP09010 

C * NSR => NUM8ER OF ANGLES PER REVOLUTION * UFP09020 

C * K => ANGLE INCREMENTS + 1 * UFP09030 

C * VX(KK) *> WHERE IME-AVERAGEO AXIAL VEL STORED * UFP09040 

C * AMP (1,1) => Z EROETH- HARMON I C AXIOAL VELOCITY * UFP09050 

C * COMPUTED BY FOUCAL * UFP09060 

C * * UFP09070 

C ********************** ***************** *************** UFP09080 

C UFP09090 

C UFP09100 

C UFP09110 

C *** COMMON TO CONTAIN GEOMETRIC PARAMETERS UFP09120 

C UFP09130 

COMMON / GEOCOM / MM . MN , NN, NW . NTORL , X( 1 2 . 1 O) . Y ( 1 2 . 10 ) , Z ( 1 2 , 10) , UFP09140 

& XTIP( 10, 1 1 ) ,YTIP( 10. 1 1 ) . ZTIP( 10, 1 1 ) .XH( 101 .7) . YH( 101 ,7) . UFP09 150 

& ZH( 101 ,7) ,XW(21 , 10) , YW(21 . 10) , ZW(21 , 10) . RH UFP09 160 

UFP09 170 

*** COMMON TO CONTAIN SINGULARITY STRENGTHS UFP09180 

UFP09 190 

COMMON / SINCOM / GT ( 10 , 10 . 60) . GB ( 10 . 9 , 60) , GTW( 10) , SB ( 10 . 9 ) , GTV UFP09200 
& , GAM (9, 60) , GMEAN( 9 ) , 0KFW(2O) , H8WX , H8WY , H8WZ UFP092 10 

UFP09220 

*** COMMON TO CONTAIN FIELO POINT DATA, NBLAOE , I DENT , NSR, AND VSR UFP09230 

UFP09240 

COMMON / FPCOM / NBLAOE , NSR , XFP , RFP , TZFP , YBAR , ZBAR , UFP09250 

& 08LA0E ,OTFP, I DENT ( 18) .TFP(60) .VSR, I PLOT, 01 AM UFP09260 

UFP09270 

*** COMMON TO CONTAIN VELOCITY OATA UFP09280 

TJFP09290 

COMMON / VELCOM / VIX ( 60 ) , VI T ( 60 ) , VI R ( 60 ) , U ( 60 , 3 ) , NSAMP UFP09300 

UFP09310 

*** COMMON TO CONTAIN HARMONICS COEFFICIENTS UFP09320 

^ UFP09330 

COMMON / HARCCM / A ( 1 5 , 3 ) , B ( 15 , 3 ) . AMP ( 1 5 , 3) , PH( 1 5 , 3 ) . NH UFP09340 

UFP09350 

01 MENS I ON VX(1) UFP09360 

UFP09370 

*** RUN THROUGH ANO FPSTEP FOR EACH TIME STEP UFP09380 

UFP09390 

00 lOO K- 1 , NSR UFP09400 

UFP094 10 

*** FPSTEP CALCULATES VELOCITY INOUCED AT A FIELO POINT BY ONE UFP09420 

*** 8LA0E AT TIME STEP K UFP09430 

UFP09440 

100 CALL FPSTEP ( K ) UFP09450 

C UFP09460 



non oooo on 



*** SUMVEL SUMS FIELD POINT VELOCITIES INDUCED BY EACH BLADE 
CALL SUMVEL 

*** FOUCAL GENERATES HARMONIC COEFFICIENTS FOR INDUCED FIELD POINT 
VELOCITY AS PROPELLER ROTATES. 

CALL FOUCAL 

*** STORE TIME-AVERAGED ( ZERO -HARMONI C ) VELOCITY IN VX 

VX ( KK ) = AMP ( 1,1) 

RETURN 

END 



UFP09470 
UFP09480 
UFP09490 
UFP09500 
UFP095 10 
UFP09520 
UFP09530 
UFP09540 
UFP09550 
UFP09560 
UFP09570 
UFP09530 
UFP09590 
UFP09600 
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SUBROUT I NE SMOOTH ( NSR . NR , VX , THETA , UA ) 

C 
C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

DIMENSION VX(80) , UA ( 60 , 9) , THETA ( 80 ) . TH( SO) . AVX(250) 

*•* COPY LAST 15 STEPS IN REV INTO FIRST 15 ENTRIES IN VX ANO THETA 

00 100 K- 1,15 

KNSR=K+NSR 

THETA ( K ) =THET A ( KNSR ) - 360 . 

100 VX ( K ) = VX ( KNSR ) 

*** COPY FIRST 5 TIME STEPS IN REV INTO 5 ENTRIES AFTER REVOLUTION 

00 200 K-16,20 

KNSR =K + NSR 

THET A ( KNSR ) =THET A ( K ) + 360. 

200 VX ( KNSR ) = VX ( K ) 

*** SET UP EVENLY SPACED ANGLES 

DO 300 K* 1 , NSR 

300 TH ( K ) = 360 . -FLOAT (K- 1 ) /FLOAT (NSR) 

901 FORMAT ( 8F 10.5) 

902 FORMAT (1615) 

FIND FIRST POSITIVE ANGLE, SAVE INDEX 



SUBROUTINE SMOOTH REAOS INOUCED VELOCITIES AT ONE 
RAO I US FROM PIEDAT. IT CALCULATES BY 
INTERPOLATION A SERIES OF NSR EVENLY SPACEO 
AXIAL INOUCED VELOCITIES SUITABLE FOR INPUT INTO 
UNSTEADY PIEWAKE. 

NSR => NUMBER OF ANGLES PER REVOLUTION 

NR => RAO I US NUMBER 1*HUB, ETC. 

VX(1 - 5) *> RAW VELOCITIES 5 STEPS BEFORE ZERO 
BLADE ANGLE. 

VX ( 6 - NSR + 5) => RAW ..VELOCITIES 'IN ONE REVOLUTION 
VX ( NSR+6 - NSR+20) => RAW VELOCITIES 15 TIME 
TIME STEPS AFTER ONE REVOLUTION 
=> ANGLE OF VELOCITY VX(K) 

=> SMOOTH VEL AT RAOIUS M. ANGLE K 
=> EVENLY SPACED ANGLES OF UA(K.M) 

=> NUMBER OF ANGLES FED TO UGLYOK 
*> K OF FIRST POS. ANGLE IN THETA ( K ) 

*> FIRST ANO LAST INOEX OF THET A ( K ) 
AND VX ( K ) SENT TO UGLYOK 
A VX => SPLINE CURVE COEFFICIENTS OF VX 



THET A ( K ) 

UA ( K , M ) 

TH(K) 

NIN 
KPOS 
KSTART , KSTOP 



UFP096 10 
UFP09620 
UFP09630 
UFP09640 
UFP09650 
UFP09660 
UFP09670 
UFP09680 
UFP09690 
UFP09700 
UFP097 10 
UFP09720 
UFP09730 
UFP09740 
UFP09750 
UFP09760 
UFP09770 
UFP09780 
UFP09790 
UFP09800- 
UFP09810 
UFP09820 
UFP09830 
UFP09840 
UFP09850 
UFP09860 
UFP09870 
UFP09880 
UFP09890 
UFP09900 
UFP099 10 
UFP09920 
UFP09930 
UFP09940 
UFP09950 
UFP09960 
UFP09970 
UFP09980 
UFP09990 
UFPIOOOO 
UFP 1 00 1 O 
UFP1OO20 
UFP 1 0030 
UFP 1 0040 
UFP 1 0050 
UFP 10060 
UFP 1 0070 
UF P 1 0080 
UFP 10090 
UFP10100 
UFP 101 10 
UFP 10 1 20 
UFP 10130 
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400 

500 



C 

c * * * 
c 



c 

c * * * 
c 



00 400 K- 1 , NSR 
KPOS^K 

IF (THETA(K) .GT.O. ) GOTO 500 

CONTINUE 

KSTART = KPOS - 1 

NSR2=NSR/2 

KST0P=KP0S+NSR2+ 1 

NIN=KST0P-KSTART+1 

FINE SPLINE COEFFICIENTS AND SMOOTH VEL FOR FIRST 180 DEGREES 

’S 

CALL UGLYDK ( NIN , 1 , 1 , THETA ( KSTART ) , VX( KSTART) ,0. ,0. , AVX) 

CALL EVALOK( NIN. NSR2. THETA (KSTART) , TH( 1 ) . UA ( 1 , NR ) .AVX) 

FIND SPLINE COEFFICIENTS AND SMOOTH VEL FOR LAST 130 OEGREES 

K$TART=KST0P-2 

kstop=kstart+nin- 1 
NSR2P=NSR2+1 

CALL UGLYDK (NIN. 1 . 1 . THETA ( KSTART ) . VX( KSTART) .0. .0. .AVX) 

CALL E VALDK ( NIN , NSR2 , THETA (KSTART ) .TH(NSR2P) . UA ( NSR2P . NR ) .AVX) 

RETURN 

END 



UFP 10140 
UFP10150 
UFP 1 0 1 60 
UFP10170 
UFP 10180 
UFP 10190 
UFP 10200 
UFP 102 10 
UFP 10220 
UFP 10230 
UFP 1 0240 
UFP 10250 
UFP 10260 
UFP 10270 
UFP 10280 
UFP 10290 
UFP 10300 
UFP 103 1 O 
UFP 10320 
UF P 1 03 30 
UFP 10340 
UF P 1 0350 
UF P 1 0360 
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** APPENDIX XI 



* * 



PIEVAXE Program Listing 



CV CN 



C THIS PROGRAM COMPUTES EFFECTIVE WAKES BY ONV METHOO 

C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

01 MENS I ON VO ( 60 ,11), RX (11), UX ( 60) , RA ( 1 1 ) , UA ( 60 , 1 1 ) .UASEG( 60 ) 
DIMENSION RE ( 1 1 ) , VR ( 60, 1 1 ) , VT( 60. 1 1 ) , AUASEG(44 ) , AUX(44 ) . VTE( 1 1 
DIMENSION THETA (60) ,VE(60. 1 1 ) . IOENT( 1 8 ) , UERX ( 60) . VRE ( 1 1 ) 

COMMON RX.UX.RA.UASEG.UERX.RE 

1 10 FORMAT ( 9 F 10.5) 

103 F ORMAT ( / ' AVERAGE VELOCITY *'.F5.3) 

113 FORMAT ( / ' VOLUMETRIC AVERAGE NOMINAL VELOCITY *'.F5.3) 

114 FORMAT ( / ' VOLUMETRIC AVERAGE EFFECTIVE VELOCITY *' f F5.3) 

115 FORMAT ( / ' EFFECTIVE BLOCKAGE *'.F7.3,' PERCENT ' ) 

111 FORMAT ( 1 8 A4 ) 

10 FORMAT ( 1615) 

OO FORMAT ( 8F 10. 5 ) 

C 

C *** REAO IN NOMINAL WAKE OATA. SAVE RAO I AL ANO TANGENTIAL 



THE ORIGINAL PIEWAKE WAS WRITTEN BY R. VAN HOUTEN 
IN ABOUT 1981. IT IS BASEO ON THE OET NORSKE 
VERITAS METHOO OF COMPUTING EFFECTIVE WAKE. THIS 
PROGRAM HAS BEEN MOOIFIEC BY R. JAMISON IN APRIL. 
1982. TO ACCOMOO AT E UNSTEAOY PROPELLER - I NDUCEO 
AXIAL VELOCITIES FROM UFPV. BASICALLY. BOTH 
PROGRAMS OIVIOE THE WAKE INTO A NUMBER OF PIE- 
SHAPEO SEGMENTS. EACH SEGMENT IS ASSUMEO TO 
CONTRACT IN AN AX I SYMMETRIC MANNER AS GIVEN BY 
T. HUANG. 



NP => 

NX *> 

RX ( M ) «> 

VO(K.M) => 

VT(K.M) => 

VR(K.M) => 

VTE.VRE »> 

UASEG(M) => 

UX ( M ) => 

AUX. AUASEG=> 
= > 
= > 
= > 
*> 
= > 
*> 
*> 
= > 
= > 
«> 



RE ( M ) 
THETA ( K ) 
NA 

RA ( M ) 
UA(K.M) 
VE(K.M) 
VL 

VLAVN 

VLAVE 

BLOCK 



NUMBER OF PIE-SHAPEO SEGMENTS 
NUMBER OF NOMINAL WAKE RAOII 
NOMINAL WAKE RAOIUS M (UP TO 11) 
NOMINAL AXIAL VELOCITY, SEG K. RAO M 
NOM TANG VEL. SEG K. RAO M 
NOM RAO VEL. SEG K, RAO M 
VT & VR INTERPOLATED AT EFF WAKE RAO 
UA IN A PIE - SEGMENT . ALSO USEO TO 
STORE A PIE -SEGMENT OF VT 
VO IN A PIE-SEGMENT. ALSO USEO TO 
STORE A PIE-SEGMENT OF VR 
SPLINE COEFF. FOR VR . VT . RESP. 

RAOII OF EFFECTIVE WAKE VELOCITIES 
ANGLE IN OEGREES OF PIE SEG K 
NUMBER OF PROP - I NO . AXIAL VEL RAOII 
PROP-INO AXIAL VEL RAOIUS M 
PROP INO AXIAL VEL. SEG K, RAO M 
EFFECTIVE WAKE VEL. SEG K. RAO M 
AVERAGE VELOCITY 

VOLUMETRIC AVERAGE NOMINAL VELOCITY 
VOLUMETRIC AVERAGE EFFECTIVE VELOCIT 
EFFECTIVE BLOCKAGE 
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PI E00010 
PIE00020 
PIE00030 
P I E 00040 
PIE00050 
PIE00060 
PI E 00070 
PIE0O080 
PIE0C090 
PIEOOIOO 
PIEOOI 10 
PIE00120 
PI E00130 
PI E00140 
PIEOOI 50 
PIE00160 
PIE00170 
PIE00180 
PIE0O190 
PIE0020C 
PIE00210 
PIE00220 
PIE00230 
PI E00240 
PIE00250 
PI E00260 
PIEOO270 
PIE00280 
PIE0O29O 
PIE00300 
PIE00310 
PIEOO320 
PIE00330 
PIE00340 
PIE00350 
PIE0C360 
PIE0C370 
PIE00380 
PIE00390 
PI E00400 
) PIE0O410 
PI E00420 
PI E0O430 
PIE00440 
PIE00450 
PIE00460 
P I E 004 70 
PIE00480 
PIE00490 
PIE00500 
PIE0C5 10 
PIE0O520 
PI E0C530 



° c> o ooo noon oor> 



VELOCITIES. THEY ARE NOT NEEDEO FDR EFFECTIVE WAKE 
CALCULATIONS BUT WILL BE ECHOED FDR WKPRDC. 

REAO(IO.m) ( IOENT(I ) , 1 = 1 , 18) 

REAO( lO. 210) NX 
READ ( 10.210) NP 
READ( 10. 200) (RX(M) ,M= 1 .NX) 

DO 20 M= 1 .NX 

REAO( 10, 200) ( (THETA(K) . VD(K.M) , VT(K.M) . VR(K.M) ) .K= 1 . NP ) 
20 CONTINUE ■ . 

READ IN PROPELLER- INDUCED VELOCITIES FROM EITHER FPUV 
OR A SEPARATE OATA FILE 



900 



22 



21 



FORMAT ( / ' DO YOU WISH TD USE STEADY OR UNSTEAOY VELOCITY PROFILE? 
& /' ENTER 1 FOR STEADY. 2 FDR UNSTEADY') 

8 WR I TE ( 6 . 900 ) 

READ ( 5 , * ) ISTED 

IF (ISTED.GT.2.DR. ISTEO.LT. 1) GOTO 8 

RE AO (8,210) NA 

RE AO ( 8 , 200 ) (RA(I ) . 1*1 .NA) 

RH=RA( 1 ) 

RPROP * 1 . 

IF(ISTEO.EO. 1) GDTD 22 

RE AO ( 8 , 200) ( (UA(KK ,M) . KK = 1 , NP ) , M= 1 , NA ) 

GOTO 92 

RE AO (8, 200) ( UA ( 1 , M ) , M* 1 , NA ) 

DO 2 1 KK = 2 , NP 
00 21 M* 1 , NA 
UA ( KK , M ) =UA ( 1 ,M) 



EXTRAPOLATE NOMINAL WAKE DATA TO- HUB 



92 WRITE ( 6 , 100) 

lOO FORMAT ( / ' 00 YOU WISH TD EXTRAPOLATE VELOCITY DATA TD THE HUB?'/ 

1 ' ENTER 2 FOR NO, 1 FOR LINEAR, O FDR CONSTANT') 

R E AO ( 5 , * ) IXTRAP 
NX 1 =NX 

IF( IXTRAP. GT. 1 )G0 TO 28 

*** LINEAR EXTRAPOLATION OF NOMINAL WAKE DATA TO HUB. 

NX=NX+ 1 
OD 25 K= 1 , NP 
00 23 M= 1 , NX 1 
MR =NX - M 
MRP 1 =MR+ 1 

VT ( K , MRP 1 ) = VT ( K , MR ) 

VR ( K . MRP 1 ) * VR ( K , MR ) 

23 VD ( K , MRP 1 ) = VO ( K , MR ) 

25 VO(K. 1 ) = VD ( K , 2 ) + 1 XTRAP* ( VO ( K , 2 ) * VO (K,3))/(RX( 1 ) -RX( 2 ) ) * ( RH-RX( 1)) 
DO 26 M* 1 , NX 1 
MR*NX-M 
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PIE00830 
PIE00840 
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P I EO 1050 
P I EO 1060 
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MRP 1 = MR + 1 

26 RX ( MRP 1 ) = RX ( MR ) 

RX ( 1 ) = RH 
28 VL =0 . 

VLAVN=0. 

VLA VE =0 . 

*•* CYCLE THROUGH EACH PIE SEGMENT 

DO 50 K«1,NP 
DO 30 M«1 , NX 
UASEG(M) =UA ( K , M ) 

30 UX ( M ) = VO ( K , M ) 

WKMOD PERFORMS THE HUANG MODIFICATION ON ONE PIE-SEGMENT 

CALL WKMOD ( NX t NX ,NA, 1 , RPROP , VOL f VOLAVN , VO LAVE ) 

• INTERPOLATE VT AND VR AT EFFECTIVE WAKE RADII 

DO 60 M S 1 , NX 
UASEG(M) =VT ( K , M ) 

60 UX(M)=VR(K,M) 

CALL UGL YDK ( NX , 1 , 1 , RX , UA SEG . 0 . , 0 . , AUASEG) 

CALL UGL YDK ( NX , 1,1, RX , UX , O . , O . , AUX ) 

CALL EVALDK ( NX , NX , RX . RE , VTE , AUASEG ) 

CALL EVALDK ( NX , NX , RX , RE , VRE , AUX ) 

DO 70 M* 1 .NX 
VT(K, M)=VTE(M) 

70 VR (K f M)=VRE(M) 

ACCUMULATE VOLUMETRIC AVERAGE DATA 

VL= VL+VOL 

VL A VN= VLA VN+ VOLAVN 
VL AVE = VLA VE+ VQLA VE 
DO 40 M- 1 . NX 
40 VE ( K , M ) = UERX ( M ) 

50 CONTINUE 

COMPUTE VOLUMETRIC AVERAGES 

VL = VL/NP 

VLA VN- VLA VN/NP 

VLA VE = VL A VE /NP 

SLOCK =( 1 . -VLAVE/VLAVN) • 100. 

*** WRITE EFFECTIVE WAKE DATA IN FILE 11 TO BE PROCESSED BY WKPROC 

WRITE ( 1 1 , 1 1 1 ) ( IDENT ( I ) , I- 1 , 18) 

WRITE( 1 1 ,210) NX 

WRITE( 1 1 ,210) (NP ,M«1 .NX) 

WRITE! 1 1 , 200 )(RE(M).M= 1 . NX ) 



P I EO 1070 
P I EO 1080 
PI EO 1090 
PIE01 100 
PIE01 1 10 
P I EO 1 120 
PIE01 130 
PIE01 140 
PI EOI 150 
PIE01 160 
PIE01 170 
PIE01 180 
PIE01 190 
PI EOI 200 
P I EO 1210 
P I EO 1 220 
PI EO 1 230 
P I EOI 240 
P I EO 1 250 
P I EO 1 260 
P I EO 1270 
P I EO 1 280 
P I EO 1 290 
P I EO 1 300 
PIE013 10 
PI E01320 
P I EO 1 330 
P I EO 1 340 
P I EOI 350 
P I EOI 360 
P I EO 1370 
P I EO 1 380 
P I EO 1 390 
P I EO 1 400 
P I EO 1410 
PI EOI 420 
P I EO 1 430 
P I EO 1440 
P I E0 1 4 50 
P I EO 1 460 
P I EO 1470 
P I EO 1 480 
P I EOI 490 
P I EO 1 500 
P I EO 1510 
P I EO 1 520 
P I EO 1 5 30 
P I EO 1540 
P I EO 1 5 50 
PIE01560 
P I EO 1 570 
PIE01580 
P I EO 1 590 
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DO 120 M 


= 1 .NX 


PI E01600 


WR I TE ( 1 1 
120 CONTINUE 


. 200) ( (THETA(K) . VE(K , M) , VT(K.M) , VR(K ,M) ) ,K= 1 , NP ) 


P I EO 1610 
PIE01620 
P I EO 1 630 


•** WRITE AVERAGES TO TERMINAL 


P I EO 1 640 
PI E01650 


WR I TE ( 6 , 


103 )VL 


PI E01660 


WR I TE ( 6 . 


1 1 3 ) VLAVN 


P I EO 1670 


WR I TE ( 6 , 


1 14 )VLAVE 


P I EO 1 680 


WR I TE ( 6 , 

STOP 

ENO 


1 15)BL0CK 


PI E01690 
PI £01700 
PI E017 10 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



10 

10 



* * * 



SUBROUTINE WKMOD ( N . NX f NA , IOW, RPROP . VOL , VOL A VN , VOL A VE ) 



SUBROUTINE WKMOD CALCULATES THE THOMAS HUANG 
CONTRACTION OF THE NOMINAL AND INDUCEO WAKES 
PASSED TO IT THROUGH COMMON. IT ALSO COMPUTES 
VOLUMETRIC AVERAGES. 



NX 

RH 

RT 

UX.UXR 
UE , UEXR 
UP 

UA , UAR 
IOW * 1 
ITER 
ITERAT 
B.C.O.F 
AUA, AUX, 

RPROP 
RE ( M ) 
VOL 
VLAVN 
VLAVE 



= > 
3 > 
*> 
= > 
= > 
= > 
= > 
= > 
= > 
= > 
= > 
AUE- 

= > 
= > 
*> 
= > 
= > 



NUMBER OF NOMINAL RADII 
RAO I US OF HUB 

RAO I US OF OUTERMOST NOMINAL VELOCITY 
NOMINAL WAKE VELOCITIES 
EFFECTIVE WAKE VELOCITIES 
APPARENT WAKE * EFF. WAKE + I NO . WAKE 
INOUCEO VELOCITIES 
NO TUNNEL CORRECTIONS 
OPEN WATER ITERATION INDEX 
HUANG CORRECTION FACTOR INDEX 
HUANG CORRECTION FACTORS. SEE PAPER 
> SPLINE COE FF. FOR INOUCEO. NOMINAL. 
AND EFFECTIVE VELOCITY FIELOS. RESP 
RADIUS OF PROPELLER 
SPECIFIED EFFECTIVE WAKE RADII 
AVERAGE VELOCITY 

VOLUMETRIC AVERAGE NOMINAL VELOCITY 
VOLUMETRIC AVERAGE EFFECTIVE VELOCITY 



RX( 1 1 ) ,R( 1 1 ) ,RA( 1 1 ) , UX ( 60) , UXR ( 60 ) , UA ( 60 ) . UP ( 60 ) 
RE ( 1 1 ) ,UAR( 60) . RP( 1 1 ) 

AUX(240).AUA(240) 

AUE(240) . UE RX ( 60 ) 



DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
COMMON RX,UX,RA,UA,UERX,RE 

USE R ( I ) FOR NOMINAL RADII 



N = NX 
N 1 = N- 1 
RH= RX ( 1 ) 

RT = RX ( NX ) 

INTERPOLATE TO FINO NOMINAL VELOCITIES AT NOMINAL RADII 

CALL UGL YOK ( NX , 1. 1,RX,UX,0. ,0. . AUX) 

DO 10 1= 1 ,N 

R( I )*(RT-RH)/(N- 1 ) * ( I ~ 1 )+RH 
R ( I ) S RX ( I ) 

CALL EV ALDK ( NX , N , RX , R , UXR , AUX ) 

INTERPOLATE TO FINO INOUCEO VELOCITIES AT NOMINAL RADII 
CALL UGLYDK ( NA , 1 . 1 , RA . UA , O . , O. .AUA) 



PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

PI 

, UE ( 60 ) PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 



E01720 
E01730 
E01740 
E01750 
E01760 
EO 1770 
E01780 
E0 1 7 90 
EO 1 800 
E018 10 
E0 1 8 20 
E0 1 8 30 
EO 1 840 
EO 1 850 
EO 1 860 
EO 1870 
EO 1 880 
EO 1 890 
EO 1 900 
EO 1910 
EO 1 920 
EO 1 930 
EO 1940 
EO 1 950 
E01960 
EO 1 970 
EO 1 980 
EO 1 990 
E 02 OOO 
E02010 
E02020 
E02030 
E02040 
E02050 
E02060 
E02070 
E02080 
E02090 
E02 lOO 
E02 1 10 
E02 120 
E02130 
E02 140 
E02 150 
E02160 
E02170 
E02 180 
E02 190 
E02200 
E022 10 
E02220 
E02230 
E02240 
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CALL EVALOK(NA,N,RA , R , UAR , AUA ) 

901 FORMAT (3F 10.5) 

-** 0UTSI0E PROPELLER DISK. NOMINAL WAKE = EFFECTIVE WAKE 

I UPD AT =0 
ITER= 1 
I TE RAT - 1 
UE ( N ) =UXR ( N ) 

20 UP(N)=UE(N)+UAR(N) 

*** ESTIMATE EFF. WAKE MARCHING FROM OUTS'lOE TO HUS 

DO 30 1= 1 , N1 
K 2 N- I 

UE ( K ) -SORT ( (UE (K+ 1 ) + (UAR(K+ 1 )+UAR(K) )/2. ) * * 2 + UXR (K)**2“UXR(K+ 1 ) * 
1 )-(UAR(K+ 1 )+UAR(K) )/2. 

UP(K)=UE(K)+UAR(K) 

30 CONTINUE 

COMPUTE HUANG'S CORRECTION FACTORS 

RP ( 1 )=RH 
VOL s 0 . 

00 40 K= 1 , N 1 
B - 2 *UP ( K+1 )+UP(K) 

C--RP(K)*( UP ( K + 1 )-UP(K) ) 

F = ( R ( K+ 1 ) *-2-R(K)**2) * (2*UXR(K+ 1 )+UXR(K) ) -R ( K ) * ( UXR (K+ 1 ) -UXR(K) ) 
1 R ( K+ 1 ) - R ( K ) ) 

VOL= VOL+F 

0 = -RP(K)**2*(UP(K+1 ) + 2 *UP ( K ) )-F 

*** FIND RADIUS FOR (K+1)ST APPARENT WAKE 

IF ( (C*C-4*B*0) .LT.O. ) GOTO 80 
40 RP ( K + 1 ) * ( -C+SQRT (C**2-4*B*0) )/2 ./B 

. IF(ABS(RP(N)-RT).LT.(.0001*RT))GO TO 60 
I UPOAT ~0 

*** I0W=0 => TUNNEL CORRECTIONS 

*1 => OPEN WATER (I.E., INFINITE FLUID) 

IF( IOW. EO. 1 ) GO TO 60 
IF ( ITER.NE. 1 ) GO TO 50 
I TER*2 

UENOLO=UE ( N ) 

RPNOLD=RP(N) 

UE(N)=UE(N)*.98 
GO TO 20 
50 CONTINUE 
ITERMTER+1 

IF( ITER. GT. 10) GO TO 80 
UEN=UE (N ) 



PIE02250 
PI E02260 
PIE02270 
PIE02280 
PI E02290 
PI E02300 
PI E023 10 
PI E02320 
PI E02330 
PIE02340 
PIE02350 
PIE02360 
PIE02370 
PIE02380 
PI E02390 
-2PI E02400 
PI E024 10 
PIE02420 
PI E02430 
PIE02440 
PIE02450 
PI E02460 
PIE02470 
PIE02480 
PIE02490 
PIE02500 
PIE02510 
r (PI E02520 
PIE02530 
PIE02540 
PIE02550 
PIE02560 
PIE02570 
PIE02580 
PIE02590 
PIE02600 
PI E026 10 
P I E02620 
PIE02630 
PIE02640 
PIE02650 
PIE02660 
PIE02670 
PIE02680 
PIE02690 
PI E02700 
PIE02710 
PIE02720 
PIE02730 
PIE02740 
PIE02750 
PIE02760 
PIE02770 



- 95 - 



ooo ooo ooo ooo o o o o ooo ooo ooo 



UE(N)=UE(N) + (RT-RP(N) ) / ( RP ( N ) - RPNDLO ) * ( UE ( N ) -UENDLD ) 
RPNDLO=RP(N) 

UENDLD=UEN 
GO TO 20 

*** ITERATE ON HUANG CORRECTION FACTORS, 3 TIMES FDR OPEN WATER 

60 IF( IUPDAT. EO. 1 )GD TO 70 

I F ( IDW . EO . 1 . AND . I T ERAT . EO . 5 ) GO TO 70 
ITER 3 1 

ITERAT-I T ERAT + 1 
I F ( ITERAT . GT . 10)G0 TO 80 
IUPOAT- 1 

FIND INOUCEO VELOCITIES AT APPARENT WAKE RADII 

CALL EVALOK(NA,N,RA,RP,UAR f AUA) 

GO TO 20 
70 CONTINUE 

*** CALCULATE EFFECTIVE WAKE RAOII 

NXM 1 = NX~ 1 
00 IOO M 3 1 , NX 

100 RE ( M ) 3 RH + ( RPRDP - RH ) * FLOAT ( M- 1 )/FLDAT ( NXM 1 ) 

IF OUTERMOST APPARENT RADIUS IS INSIDE PROPELLER OISK, THEN 
LINEARLY EXTRAPOLATE TO PROPELLER RADIUS 

IF (RP(N) .GE. RPRDP) GOTO 110 

UE ( N ) =UE ( N 1 ) + ( UE ( N ) -UE(N1 ) ) *( RPRDP -R(N1 ) ) / ( RP ( N ) - RP ( N 1 ) ) 
RP(N) 3 RPRDP 

*** INTERPOLATE EFFECTIVE WAKE VELOCITIES AT SPECIFIED RADII 

110 CALL UGLY OK (N. 1, 1 ,RP,UE,0. ,0. ,AUE) 

CALL EVALDK ( N , NX , RP , RE . UERX , AUE ) 

*** COMPUTE VOL = AVERAGE VELOCITY 

VOL 3 VOL/ 3./(RT**2- RH* * 2 ) 

*** COMPUTE VOLUMETRIC AVERAGE NOMINAL VELOCITY 

CALL INTEOK ( NX , RX , RH , RPROP , YOX . XYDX . XXYDX , AUX ) 

VOL AVN-XYOX *2/(RPROP**2 -RH* *2 ) 

*** COMPUTE VOLUMETRIC AVERAGE EFFECTIVE VELOCITY 

CALL I NT EOK ( N , RP , RH , RPROP . YOX . XYDX , XXYDX , AUE ) 

VDLAVE*XYDX*2 . / ( RPROP * * 2 -RH* -2 ) 

RETURN 
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TUNNEL CORRECTION OR HUANG CORRECTION ITERATIONS FAIL TO CONVERGE 

STOP 

END 
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